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Chapter 1 
mTBODDCTlOh 


1 . I Motivation 

The turbulent shear layer is an important flow structure in most 
industrial combustors. Premixed flames are stabilized by heat 
recirculation involving the recycling of hot products %rhich are mixed 
with cold v^actants in a turbulent shear layer. Diffusion flames are 
controlled by the rate at which fuel and oxidizer mix in a turbulent 
shear layer. 

There has been considerable interest in premixed combustion 
related to gas turbine engines because of the pollutant advantages 
(Jones et el» 1978). Premixed combustion in model gas turbine 
combustors can give a twentyfold decrease in nitric oxide emissions 
(Lefebvre, 1977). A particularly simple experimental configuration 
incorporating the aspects of premixed combustion stabilized in a 
turbulent shear layer is the rearward- facing step combustor (Figure 
1.1). A premixed flame forms in a two-dimensional turbulent shear 
layer formed at the edge of a rearward- facing step. Stabilization is 
achieved by recirculation of hot products behind the step %ihich mix 
with the incoming premixed reactants in the shear layer. 

The two-dimensional, rearward-facing step combustor facilitates 
easy access for experimental viewing by advanced laser diagnostics and 
flow visualization. Measurements in turbulent, reacting flows are 
limited and necessary for comparison to turbulent models. The 
rearward-f ac ing or backward-f ac ing step flow has been the subject of 
considerable modeling efforts (Marvin, 1977, Ashurst, 1979). Also 
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recently Qtoniem et el (1980) used a rsndooi vortex technique (Chorin, 

1973) to model the rearward- facing step combustor, the 
t vKi-dimonsional ity of the rearward- facing step flow make it more 
tractable for modeling. 

In earlier work» Ganji and Sawyer (1980) described the flow field 
of the rearward- facing step combustor with high speed schlieren 
photography. They found the flow field is dominated by large scale 
structures with an ordered pattern previously unreported in combustion 
flows. The react ing, turbulent shear layer structure is similar to 
that found in turbulent mixing layers where the importance of large 
scale structures has already been established (Brown and Roshko, 

1974) . The purpose of this work is to quantify the effect of 
combustion on the shear layer dynamics. The turbulent velocity is 
mapped in detail by laser Doppler velocimetry (LDV) for both 
isothermal and reacting flows. The effect of ccxnbustion on the 
rearward- fac ing step flow properties such as the shear layer growth 
rate» entrainment rate, turbulence development, and large scale 
structure development are determined. 

I . 2 Background 

The turbulent structure of the rearward- fac ing step is largely 
preserved under combustion. A discussion of its flow field will lead 
to a better understanding of the reacting flow. The backward facing 
step flow can be separated into two regions: i) the mixing layer 
region and ii) the relaxation region (Figure 1.2). The mixing layer 
region includes the flow from the initial boundary layer at separation 
to the reattachment point. The relaxation region covers from flow 


reattachment to the point of full recovery of a turbulent boundary 
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layer. The first region noraally extends 5-7 step heights 
doimstream of the step (Eaton and Johnston, ]980a). The reattactment 
region covers fr<n& reattachment to a distance of about 30 times the 
mixing layer thickness at reattachment (Bradshaw and Wong, 1972)* 
This study focused on the mixing layer region. The mixing layer in a 
rearward- facing step flow is often called a reattaching shear layer. 
The reattaching shear layer is initially planar but then curves toward 
the wall and eventually impinges on the wall at reattachment. The 
turbulent structure of planar mixing layers is very similar to the 
reattaching shear layer. Since planer mixing layers have been studied 
extensively, their current state of understanding will be discussed 
next . 

Mixing layers . The current understanding of the nature of 
turbulent mixing layers stems from the work of Brown and Roshko 
(1974). In a series of shadowgraph photog*^aphs they showed that the 
layer is composed of large scale, two-dimensional rollup vortices. 
The vortices or eddies evolve from a Kelvin-Helmholtz instability of 
the initial shear layer (Sherman, 1976). These two-dimensional 
large-scale structures grow by entrainment of fluid and pairing as 
they are converted downstream (Winant and Browand, 1974). Pairing of 
eddies is the uniting of two neighboring eddies into one with twice 
the original size. The pairing process is a result of the inherent 
instability of two parallel line vn»-tices (Moore and Saffman, 1975). 
The growth of the eddies determines the growth of the shear layer as a 
whole . 

The process of f.rst bringing the irrotational fluid into the 
shear layer and then later ir.ixing it microscopically is given the 
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names entrainment and mixing . The two processes are found to be almost 
separate and distinct in mixing layers (Dimotakis and Bro%m, 1976). 
The large scale structure growth which is responsible for entrainment 
is mainly a two-dimensional process and not a strong function of 
Reynolds number. 

In the initial shear layer » mixing is a molecular diffusion 
process within the large scale structures. At Reynolds numbers 
greater than 10*^ based on the shear layer thickness, Konrad (1976) 
found a three-dimensional instability appears which increases mixing 
by 25%. More recently, Roshko and Bernal (1981) reported an intricate 
three-dimensional structure to the two-dimensional eddies which 
certainly enhances mixing. Thus the mixing is found to be Reynolds 
number dependent . 

The initial conditions can have a pronounced influence on the 
large scale structure formation and the subsequent shear layer growth 
and turbulence development. In a comprehensive review of mixing 
layers, Birch (1980) suggests that the shear layer flow will not 
become completely independent of the initial conditions until 
> 2 X 10^. At lower Reynolds numbers, the state of the boundary 
layer will inflverce shear layer growth rate, virtual origin, and 
turbulence development. The relationship between the initial 
conditions and the shear layer properties is still unresolved. Some 
researchers report an increase of the growth rate when the inlet 
boundary layer is tripped (Wygnanski and Fielder, 1970, Batt, 1975) 
while others report a decrease (Browand and Latigo, 1978). Husain and 
Hussain ( 1979 ) tripped the boundary layer in a half jet mixing layer 
and found t..e growth rate is lower at first and higher later. 
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High free stream turbulence influences the shear layer 
development. Birch (1980) suggests that the free stream turbulence 
primarily influences the shear layer through the inlet boundary layer 
and can be ignored for levels less than 0.6Z. Chandrsuda et el (1978) 
found that at high free stream turbulence levels the two-dimensional 
eddies largely disappear and suggest that the asymptotic state of the 
turbulent mixing layer is more three-dimensional. 

Rearvai d- f ac ing step : The mixing layer region . The mixing layer 
region of the rearward- facing step flow, often called the reattaching 
shear layer, differs from the planar mixing layer. The reattaching 
shear layer is initially planar, then curves toward the wall and 
impinges on it at reattachment. This curvature is thought to be 
responsible for the rapid decrease in shear stress near reattachment 
(Eaton and Johnston, 1980a, Castro and Bradshaw, 1976). 

The recirculation zone behind the step results in variable, 
non-zero velocities that cause the velocity difference across the 
layer to vary with downstream distance. The velocity difference in a 
planar mixing layer is constant. Also the turbulence levels in the 
recirculation zone are higher (Kim, Kline and Johnston, 1978). 
Pressure driven, stress bearing fluid is swept into the recirculation 
zone from the reattachment region (Bradshaw and Wong, 1972). This 
inceases the turbulence levels in the mixing layer (Eaton and 
Johnston, 1980a). Again the initial conditions influence the flow 
development. Besides effecting the mixing layer properties discussed 
earlier, variation of the boundary layer state changes the 
reattachment length (Eaton and Johnston, 1980a). A comprehensive 


review of backward- fac ing step research is given by Eaton and Johnston 
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Rearward- facing step combustor . The only previous experiments on 
the rearward- facing step combustor were performed by Ganji and Sawyer 
(1980). Using high speed schlieren photography they discovered that 
the reacting flow is dominated by large scale coherent structures 
found in mixing layers. Combustion is primarily confined to the 
eddies as they entrain premixed reactants and hot products. The 
spread of combustion is linked to the growth and development of the 
large scale structures. Ganji and Sawyer also described the pollutant 
and stability characteristics of the flow. They made time average 
measurements of the major gas species throughout the flow field and 
obser'^ed flashback and blowoff with high speed schlieren photography. 

1 .3 Goals 

The primary aim of the present work is to assess the effect of 
combustion on the turbulent structure of the rearward- facing step 
flow. This includes the effect of combustion on the shear layer 
growth 9 turbulence levels, reattachment length, and large scale 
structure development. To carry out this goal the following objectives 
are set down: 

1. To construct and develop a frequency shifted, laser Doppler 
velocimeter (LDV) and make unambiguous, time-resolved, velocity 
measurements in the turbulent, reacting, recirculating flow. 

2. To develop the computer hardware and software to acquire the 
instantious velocity measurement and calculate the appropiate 
turbulent quantities . 

To obtain a detailed mapping of the mean velocity and turbulent 
intensity in both the reacting and non- reacting flow. To 


3 . 
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evaluate the shear layer growth rate , entrainment rate , and 
recircluation length in both cases. 

4. To record the velocity probability density functions (PDFs) and 
determine the effect of combustion on the higher order moments. 

5. To compare the large structure formation and development in the 
reacting and non-reacting shear layers by frequency analysis of 
the LDV signal and visualization of the flow field by high speed 
schlieren movies. 
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Chapter 2 

BXPERIMERTAL APPARATUS AHD mSTRDMEHTATION 

2 . 1 The Two-Dimensional Combustor 

The two-dimensional combustor apparatus was basically unchanged 
from that described by Ganji (1979) or Ganji and Sa%iyer (198C). The 
major improvements included accurate pressure regulation, lowered 
inlet turbulence levels, and elimination of the air cooling in the 
test section. 

A schematic ol the overall system is viewed in Figure 2.1 and 
photographs are shown in Figures 2.4 and 2.10. The coordinate system 
used in the point measurements is given in Figure 2.3. A compressor 
supplied air at a maximum gage pressure of 630 kPa and a maximiim 
flowrate of 0.3 kg/s. A Fisher series 99 pilot operated pressure 
regulator maintained the air pressure within of the set value. 
The air was dried below 20Z humidity by a desiccant drier and filtered 
by a Balston air filter (Balston Filter Products, Lexington, Mass.). 
The Balston type A filters have a 99.99991 retention efficiency for 
0.6 micron particles. The air stream was split and metered by t%^ 
son^c nozzles which were calibrated at the Naval Air Rework Facility 
in Alameda, California to an uncertainty of less than .232. During 
the laser Doppler velocimetry (LDV) measurements, additional air 
containing the particle seed was added directly downstream of the 
nozzles (See Section 2.3.6). The air was recombined at a manifold and 
passed through three parallel Venturi tubes. Fuel was injected at the 
throat of each Venturi tube. 

The propane fuel was ^,upp1 ied by two 45 kg bottles of liquid 
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propane placed in atean heated water batha to aaintain a conatant fuel 
prceaurr. The fuel waa mtered by a Pisher-'Porter (No. 
FP-3/4-‘27 -k;- 10/80) or a Natheaon 30$ rotaaMter depending on the 
flowrate requireaenta. The fuel flow waa controlled by a needle valve 
do«matreaa of the rotaaetera. The fuel then paaaed thru a aolenoid 
ahutoff valve and flane arreator before aixing with the air in the 
venturi aection. 

The croaa aection of the coabuator ia aho%m in Figure 2.2. After 
the fuel and air were aixed in the venturia, they entered a one aeter 
long pnnnix aection 31 an high and 173 am wide. The aection waa 
packed with fine grade atainleaa ateel wool to reduce the inlet 
turbulence level and prevent flaahbark. A atainleaa ateel acreen (80 
aexhl waa atretched acroaa the end of the preaix aection to contain 
the packing. A braas foil aafety port waa deaigned to blowout at 
exceaaive prcaaurea. The flow converged over the backaide of a 
pcv-tfiled alep (Ganji, 1979) with a 2:’ area ratio. The 
t wo-dimenaiv^nal ahear layer formed at the edge of the 23 aan high atep. 

The teat aection waa 2^0 nm long, 173 mm wide, and 31 mm high 
givi. ’ an aapect rat iv> (channel width to atep height) of 6.9. Quartc 
windows (12.7 mm thick. Corning 79AO fuaed ailica) were inatalled on 
the aides of the teat aection expoaing 13 mm of the atep and extending 
the entire length of the teat aection. Acceaa porta (4.8 omi dia.) 
were provided on the centeiline of the top and bottom platea of the 
test section at 10 mm intervals downatream of the atep. The hot 
exh.iuat gases paaaed out of the teat aection into a 100 mm diameter 
exhaust pipe. Cooling water flowed thru the rear of the teat aection 
and was spr.iyed into the exhatist . A disc in the exhauat line waa 
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rotated to adjuat the test section pressure uhich was set at one 
atnosphere for all tests. 

A safety system was designed to shut down the combustor in case 
of hasardous conditions. Detectors were placed throughout the system 
(see Figure 2.1) and shut the combustor down under the following 
circumstances: 

Detector Harardous condition 

1. Flame detector Flame in premix section (flashback) 

2. Over pressure Excessive pressure in premix section 

3. Water pressure Cooling water loss 

4. Under temperature Exhaust temperature drop (flameout) 

5. Air .:iow Air flow loss 

6. Fuel leak Propane leak in room 

The steam heated fuel baths were equipped with a independent safety 
system that turned off the steam and rang an alarm if the water bath 
temperature became too high. 

To ignite the combustor » a high voltage ignitor rod was inserted 
into the shear layer from the top plate. The ignitor was located 30 

mm downstream from the step and was normally withdrawn from the test 

section after ignition. 

Chromel alumel thermocouples were installed throughout the system 
and connected to a single selector switch. A digital thermometer 
(Omega model 2160A-K) displayed the temperature with an accuracy of 
^IK. Thermocouples measured the gas temperature at the sonic noaxles^ 
fuel rotameters, and the test section entrance. Pressure transducers 

(Senso Metri'' model SP91) sensed the pressure upstream of each norzle 

and fuel rotameter with an accuracy of It of full scale. A single 
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digital panel meter indicated the pressure readings* 

2.2 Data Acquisition System 

A PDP 11/34 computer system was utilized to calculate and store 
the overall flow conditions of the combustor luring the experiments. 
The software package contained a... the calibratii^n curves of the fuel 
rotameters, sonic nozzles, pressure transducers, and the air rotameter 
(only used during LDV measurements). All the pertinent pressures and 
temperatures were directly keyed into a CRT terminal. The desired 
equivalence ratio, ^ ,and reference flow velocity, Uo, (average inlet 
velocity) were input and the program calculated the corresponding 
settings for the sonic nozzles and fuel rotaaeter. A sismary of the 
combustor data stored by the computer program is shown in Table 2.1. 
The system temperatures and pressures were constantly updated to 
maintain the desired equivalence ratio and inlet flow velocity. ^e 
accuracy of the inlet velocity and the equivalence ratio settings %«ere 
1.3S and 6Z respectively. 

2 . 3 Schlieren System 

An important aspect of the experiments was a series of high speed 
schlieren movies and long exposure photographs taken of the turbulent 
free shear layer. A conventional "z" configuration was used in the 
experiment (Figure 2.3). The output of a 1000 watt xenon arc lamp 
(Oriel Optics, model C-60-50) was collimated by a 3.94 M focal length, 
0.3 M diameter spherical mirror. The collimated light passed through 
the test section and was refocused by a second identical mirror onto 
the schlieren stop plane. A lense imaged light emanating from the test 
section onto the camera film plane. A 102 x 127 mm box camera was used 
for th«* long exposure (1/30 sec) photographs. The high speed movies 
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%rere recorded by a Hycw caera (Redlake Corp., aodel 41-0004) at 
rates ranging fron 3000-8000 fraoies/sec. 

':k>vie8 visualizing the entire flow field and enlarged high 
resc*ution movies of just the initial portion of the flow field were 
ott .Tied. Complete specifications of the schlieren optics are given 
in Table 2.2 and Figure 2.3. 

\ typical polar diagra of the deflected light at the schlieren 
stop plane is shown in Figure 2.6. The distinct vertical and 
horizontal deflections result from taperature gradients on the 
v:ombu;tor surfaces. Thermal boundary layers on the top of the step, 
and the top and bottom plates of the combustor produced the bright 
vtrtical lines. Light deflected from thermal gradients near the front 
of the step generated the horizontal line. The broad, disperse 
pa tern on the photograph resulted from index of refraction gradients 
in the shear la' . In order to visualize the shear layer, the 
schlieren stops were taylored to highlight the disperse deflections 
from the shear layer and minimize interference from the wall boundary 
layers . 

The schlierin photography was only used for visualization of the 
shear layer. o attempt was made to deduce temperature information 
from th< photographs. Ganji and Sawyer (1980) give a detailed 
disct sion of light refraction in a combusting media. 

? . Hot <Jire Anemoirelry 

Although roost of the velocity measurements were carried out using 
laser Do pier velocimetry, the superior frequency response and ease of 
op'. at ion of hot wire anemometry made the technique useful for the 
characterization of the entry flow. The hot wire measurements were 
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limited to the entry region uhere the flow direction was unoftbigjous. 
Hot wire measurements in the entry region also served as a check on 
ihe LDV velocity measurements. 

The hot wire equipement included a TSl model 1010 Heat Flux 
System and a model 1005B Linearizer. Wollaston wire (lOZ rhoduim» 90Z 
platinum) was soldered across a TSI 1210 **u** probe. The sensor had an 
active length of 1 mm and a diameter of 8 microns. The hot wire probe 
was calibrated in the low speed wind tunnel facility. A typical 
calibration curve is shown in Figure 2.7. 

Both mean velocity and rms turbulence measurements were made in 
the entry flow. The hot wire voltages %iere displayed on a digital 
volt meter capable of measuring the true rms voltage. The single 
sample uncertainty (Kline and McClintock^ 1933) of the mean velocity 
measurement at a typical velicity of 10 m/s was 1.5Z. The error in 
the rms velocity was much more difficult to assess. The frequency 
response of the probe was estimated with a square wave test to he 3 
kHz. Individual LDV particles follow the flow more closely than this 
and the nns turbulence levels measured by the hot wire %iere lower than 
the LDV measurements. Also at velocities above 10 m/s the probe 
vibrated at 9.3 kHz producing specious rms turbulence of the same 
order as the actual flow turbulence. 

The frequency spectrum of the hot wire fluctuations were recorded 
by a digital fast Fourier transform analyzer (GenRad model 2312). 
Shear layers are sensitive to flow instabilities which include 
account ical resonances in the combustion tunnel. The hot wire 
spectrum le the entry flov was checked to confirm that all such 
resonances were removed (See Section 3,1 for more details). 
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2.5 Laser Doppler Velociaetry 

Due to the hostile nature of the flov environaent^ laser Doppler 
velocimetry was selected to quantify the velocity flow field. The hot 
recirculating, turbulent flow precluded the use of any other 
technique. Hot wire anemometry, widely used in turbulent shear flow 
research, cannot withstand the high temperatures or resolve the flow 
direction. Laser Doppler velocimetry (LDV)» also called laser Doppler 
anemometry (LDA), is able to obtain measurements in combusting flows 
and resolve the velocity direction while maintaining high frequency 
response and spacial resolution. LDV has the additional advantages of 
producing a linear output and causing no flow disturbance. 

The technique was first demonstrated by Yeh and Cummins (1964) 
and since that time has widely used for measurements in laminar and 
turbulent flows. Applications to non-reacting flow phenomena are 
extensive and have been revie%red by Durst et el (1972). Measurements 
in combusting systems are not as numerous and subject to the unique 
problems of high temperature environments. Index of refraction 
inhomogeneities caused by chemical concentration and temperature 
gradients lead to difficulties in achieving a stable location of the 
laser teams in the flow. Scattering particles are required for LDV 
and these must be made of refractory materials able to withstand the 
high temperatures. The special problems of design and application of 
LDV systems for combustion have been discussed at length by Self and 
Whitelaw (1976) and Chigier (1977). 

2.5.1 Principles of Operation . The LDV measurement technique is 
based on the well known Doppler effect involving the frequency shift 
of light waves scattered by a moving object. With the development of 
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lasers which provided a coherent, monochromatic source of light, this 
effect was exploited to make precision point measurements of 
velocity. Particles are required to give sufficient scattered light 
intensity and the technique is limited by the ability of the particles 
to follow the flow field. There are numerous types of LDV systems and 
these are described in various reviews and introductions to the topic 
such as Drain (1980), Durst et el (1972, 1976), Goldstein and Rreid 
(1976), Trolinger (1974), Stevenson (1977), and Whitelaw (1975). 

The two most predominate optical configurations for LDV are the 
dual beam, real fringe, differential scattering system and the single 
scatter, reference beam system (Durst et el, 1972). Due to the ease 
of alignment and better signal quality, the dual beam, real fringe 
system was chosen for this experiment. In gas flows, idiere the 
particle density is low and the number of particles in the focal 
volume are few. Drain (1972) has shown theoretically that the dual 
beam system will give a superior signal to noise ratio (SNR) than the 
reference beam system. 

A typical dual beam, real fringe optical configration is 
illustrated in Figure 2.8. The laser output is split into two beams 
and focused into the flow. An ellipsoidal volume of planar fringes is 
formed at the intersection of the two beams with a spacing, 

df * X/i 2sin(0/2)I 2.1 

where \ is the laser wavelength and 0 is the beam crossing angle. The 
beam crossing angle is given by the focal length of the focusing lens, 
f, and the beam separation, djj, (0 * sin*'d 3 /f). Particles passing 
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through the fringe pattern will scatter light with a frequency, 

Vjj e dfU 2,2 

where U is the velocity of the particle perpendicular to the planar 
fringes. The focal volume in an ellipsoid (Adrian and Goldstein, 
1971) with the dimensions (Figure 2.8): 

dx " d/cos(0/2) 

dy ■ d 2*3 

dz = d/sin(0/2) 

The waist diameter, d, is the diffraction limited size of the laser 
beam, 


d = 2Af/(iTD) 2.4 

where D is the laser beam diameter (defined at the 1/e intensity 
points) at the focusing lens (Goldstein and Kreid, 1976). 

Light scattered by particles moving through the focal volume is 
modulated at the Doppler frequency given by equation 2.2. The 
collecting lens images light emanating from the focal volume on to a 
photomultiplier tube (PMT). Signal processing deduces the Doppler 
frequency from the PMT output. 

The optical setup shown in Figure 2.8 has one drawback. A 180 
degree shift in the velocity direction will result in an identical 
Doppler frequency. This is refered to as "directional ambiguity". 
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The difficulty is overcome by shifting one of the beams to a slightly 
different i^requency %diich causes the fringes to move at a constant 
velocity (Durst and Zar^^ 1974). A particle that is stationary in the 
focal volume will scatter light at the shifted frequency, Vq* As long 
as the fringe velocity is higher than any particle velocity in the 
same direction, all velocities will result in a unique, nonzero 
Doppler shift given by. 


^D “ ^f ^ 

In this experiment an acousto-optic modulator, a Bragg cell, produced 
the frequency shift (Durao and Whitelaw 1975). 

2.5>2 The Optical System . The actual optical system designed for 
the experiment is illustrated in Figure 2.9. The output of an 
argon-ion laser tuned to the 514.5 nm line (0.3 watts) passes through 
two lenses (11 and 12) and a polarization rotator (1/4 wave plate, 
Oriel model 2566). A path compensated beam splitter (OEI, Karlsruhe, 
West Germany) outputs two beams with a 50 mm separation. The beam 
splitter is polarization sensitive and the 1/4 wave plate is rotated 
to equalize the intensity of the output beams. Also the fringe 
pattern intensity is dependent on the polarization of the beams. 
Maximum interference is obtained when the two beams are linearly 
polarized perpendicular to the plane containing them (Stevenson, 
1977). 

One beam passes through the Bragg cell and is down shifted by 40 
MHz causing the fringes to move in the opposite direction to the main 
flow. The Bragg cell (TSI model 982) is tilted to obtain proper phase 
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matching (Yarif, 1975) of the input and output beams resulting in a 
Bragg cell efficiency of 60Z. The other beam passes through a path 
compensator and a neutral density filter (60Z) to maintain equal path 
lengths and beam intensities « The laser beams are focused into the 
test section by an 82 mm diameter, 300 mm focal length lens (13). 
Optical access is provided by 12 n«.i thick quartz windows (Corning 
7940). 

The probe volume dimensions are given in Figure 2.9. The probe 
dimensions, position, and fringe parallelism are all subject to the 
propagation properties of Gaussian laser beams predicted by scalar 
diffraction theory (Hanson, 1973, Durst and Stevenson, 1977). For 
proper alignment, the laser cavity, lens positions, focal lengths, 
optical path lengths, and access windows must all be taken into 
account. The argon-ion laser operated in a TEMq^ mode producing a 
Gaussian shaped output beam. Lenses 11 and 12 were selected to 
position the laser waist at the beam intersection with a diameter of 
96 microns. 

In line, forward scatter collection optics were used to obtain 
the maximum SNR (Durst and Whitelaw, 1971). Particle scattering is 
governed by Mie scattering theory (particle diameter =? \ ) where the 
maximum scattering intensity is in the forward direction. The laser 
beams were blocked by beam dumps and the scattered light from the 
focal volume was collected by an f/6 lens (14) and focused on the 0.5 
mm diameter pinhole of the PMT housing. An adjustable aperture stop 
was positioned after the collection lens to minimize scattering from 
tho quart:, windows and adjust light levels on Ihe PMT to prevent 
s at ur .It 1 *>n . A 1 nm bandpass intorferonco filter was installed in the 
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light tight PMT housing to sliminate room light and flsBS aaission. 
The Pacific Photonetric aodel 3150 PMT housing incorporated «u-aatal 
tube encasement and a nickle coated exterior to shield sgainst RF 
noise. An EMI 9^blA/B PMT detected and amplified the scattered light 
and output the Doppler signal for processing. 

The laser, transmission and collection optics, and the PMT 
housing were all secured to a single optical bench. The b*nch top was 
mounted to a milling table %ihich allowed the test section to ’>e 
scanned in all three directions. The LDV optics, laser, and milling 
table are shown in a photograph of the experimental apparatus in 
Figure 2.10. Most velocit> profiles were taken vertically across the 
shear layer and a stepping motor drive was installed to easily move 
the table in the y direction. A dial indicator with a resolution of 
0.01 mm displayed the y position. The y direction contained the steep 
velocity gradients and the milling table could be positioned to J^.OS 
mm. The and z directions had a positioning accuracy of ^0.5 mm. 

2.5.3 Signal Processing . The signal processing system for the LDV 
measurements is shown in Figures 2.10 and 2.11. The output from the 
PMT is amplified by a 140 MHz bandpass, xlOO gain preamplifier 
(Pacific Photometric model AD-6). The frequency shift (40 MHz) of the 
single Bragg cell is undesirably large. Therefore the shift is 
effectively reduced by the ''do%nunix" circuit of the TSI model 985 
driver/dovnshifter . A 40 MHz oscillator drives both the Bragg cell 
power amplifier and the downmix circuit. The effective frequency 
shift is selectable from 10 kHz to 20 Mhz in a 1, 2, 5, sequence. The 
frequency shift had a calibration accuracy of 0.00252. 

A Macrodyne model 2096 was selected to analyse the Doppler 
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frequency. Period counters are prefered over frequency trackers for 
analysis of single rarticle U>V signals found in turbulent gas flows. 
Under low duty ratios, they can be operated with wide input 
bandwidths allowing a wider range of particle velocities than 
frequency trackers (Humphrey et el 1975). 

The block diagrams for the Hacrodyne 2096 moduals are given in 
Figure 2.12 and 2.13. In the Front End Detector (FED), the Doppler 
signal is first amplified (xl or xlO) and filtered. The selectable 
low pass and high pass filters (0.5, 2, 4, 8, 16, and 32 MBs, 40 
db/decade) serve two functions, removing the pedestal and filtering 
out background noise. A threshold level descriminator on the filtered 
output prevented processing of residual noise. An overload level 
discriminator on the amplifier output eliminated large particle or 
multi-particle bursts. 

The burst signal enters a zero crossing detector coupled with a 
Schmidt trigger (Fig ire 2.12). The Schmidt trigger outputs a pulse 
everytime ‘be signal crosses zero and behaves like a sinewave mth an 
amplitude exceeding the threshold level. 

The Logic and Output Display (LOG) nodule performs the 5/8 
validation test (Figure 2.13). The time for 5 pulses and the time for 
8 pulses is compared. If they are not within a preset accuracy, the 
signal is not accepted. The Kacrodyne has either a 5/8 or 10/16 
comparator for the validation test. The 5/8 comparator at an accuracy 
setting of 2 was used in the experiment. This required that the two 
periods compare to about IZ accuracy. If the burst was validated a 
"sync" pulse was output. A validated burst is often described as a 
"realization". The sync pulses were externally counted to display the 
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instantaoeous data rate (Figure 2.11). The period for 8 pulaea (burst 
period) was tined by a 10 bit (1024 count) counter with variable 
ranges to cover 1.2 kHz to 100 MHz input frequencies. 

The counter had a resolution of 2 ns. At a typical Doppler 
frequency of 8 MHz using the 5^8 comparator, the burst ^ ->d 's 
measured to an accuracy of 0.2Z. The 10 bit counter outputs toe 
period with the same accuracy. A typical burst period of 500 counts 
would have an a«.curacy of 0.2Z. 

A multiplexer and external clock (0.1 or 1 MHz) were added to 
output the time between realizations or the data period (Figure 
2.13). The burst period, counter range, comparator setting, 
comparator accuracy, and data period were multiplexed over a 16 
parallel line interface to a PDF 11/34 computer. An internal D/A in 
the Macrodyne was connected to a digital volt meter (DVM) for an 
continuous display of the burst period (Figure 2.11). 

2.5.4 Velocity Data Acquisitio n and Reduction . A typical LDV data 
record is shown in Table 2.3. The LDV system parameters and typical 
values are indicated. In addition to the data directly accessed by 
the parallel line interface, information on the record length, fringe 
spacing, Bragg cell shift, and probe location were keyed into the 
computer. The LDV computer control program collected a 2048 point 
record on comand from the computer terminal. The instantaneous 
velocities were calculated accoi *ing to. 


“1 • t'^Di " Vo)/df 


%rhere is the Doppler frequency of an individual burst. Two types 
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of velocity everesea were calculated, nwerical and tiaM integrated 
averagea. The otaMrical averagea are given by: 


S “ 


2.7 


.TBS • I 1/N|(U1 -U)*!** 2*® 

%fhcre M i« the record lengthy QT is the mean velocity and Uxme 
rms velocity. These averages vhich are calculated from individual 
particle bursts are subject to velocity biasing discussed by 
McLaughlin and Tiederman (1973). In turbulent flov the higher 
velocity particles have a higher probability of crossing the focal 
volume than the lowr velocity particles resulting in a bias. 
Dimotakis (1976) suggested that in high data rate regimes the biasing 
error could be circumvented by the use of the time integrals, 

TI, - l/T/u(t) dt 2.9 


. , A 2- 10 
«rn.«t - ll/T/ (ui-ut)dtl^ 

where and sre the tine integrated mean and ms velocities, 
respectively. The integrals were in.*grated nunerically by the 
rectangular rule and tine weighted averages of the data r^rord were 
stored at all neasurenent locati,jns. 


The probability denaity function (PDF 


ty wan also 
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calculated and atored. The velocity tecord waa divided into 72 equal 
intervals (bine) between the •axiaiB and suniniiai velocitiea. The 
nunber of velocity points in each bin «ao computed to deteraine the 
PDF. 

Again a time correction was applied to avert the velocity biasing 
in the PDF at high data rates. Each data point was weighted by the 
data period (dt). Typical PDFs in the shear layer are shown in Figure 
3.28. 

Another biasing problem occurred with the PDFs. The Macrodyne 
output is inversely proportional to the velocity, that is the hurst 
count register is in equal units of time (burst counts) not 
frequency. Equal intervals in velocity in the 72 bin PDF will not 
correspond to equal intervals in burst counts. For some high velocity 
(low burst count) intervals there will be no burst count value 
corresponding to that velocity interval. The PDF intervals must be 
weighted by the ratio of dBe/dBa idiere dBe is the exact width of the 
velocity interval in burst counts (ie 2.51) snd dBa is the possible 
ntmber of burst count values that fall in that interval (ie 2). All 
the PDFs were weighted in this fashi.on. This did not completely 
eliminate the biasing problem. The correct method is to store the PDF 
of the burst period and produce velocity PDFs with equal intervals of 
burst period and not velocity. 

Spectral analysis of the Doppler signal was also performed. 
Large data records (50,000 - 200,000 points) involving sample periods 
of a minute or more were required to obtain frequency spectra with 
sufficient signal quality. The PDP 11/34 had no provision for storage 
of such large records and the on line fast Fourier transform (FFT) 
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software ms very slow. Bence a FFT analyser (GenRad aodel 2512) 
processed the analog output fro* the burst period counter. Up to 512 
FFTs %iere averaged by the analyser to obtain the spectrm. A typical 
spectruB is shown in Figure 3.46. This is the spectriss of the a.c. 
coaiponent of the burst period rather than the Doppler burst 
frequency. For low turbulence levels the relationship is linear, but 
for high turbulence levels the anplitude of the spectra will be 
biased. The aajor use of Che spectra was to quantify the frequency of 
the large scale structures in the shear layer and the value c.f the 
spectral density was unimportant. 

2.5.5 Data Verificat ion . The LDV system had to be checked to 
ascertain whether the system was generating the correct velocities and 
turbulence levels. The greatest source of systematic error involved 
adjustment of the Macrodyne front panel settings: the threshold and 
overload discriminators, the lotrpass and highpass filters, and the 
counter range. 

The threshold setting was determined by Che amount of noise 
measured at the "Filtered Output" (Figure 2.12). The threshold level 
was set such that the noise would not be validated. First, the 
threshold level was calibrated (Figure 2.14). The minimun validated 
sinewave anplitude was measured as a function of the threshold level. 
For a given threshold setting, if the Doppler burst is above that 
amplitude and fulfills the other Macrodyne validation criterion, it 
will be validated. The threshold level could simply be set above the 
noise level given by this curve. 

This was found to be too strict of a criterion. The turbulence 
intensity and data rate as a function of the threshold setting is 
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•hum in Figure 2.15* At lov threshold settings the noise v^uld be 
vslidated end the turbulence intensity uould depend strongly on the 
threshold setting. As the threshold level mas increased the 
turbulence intensity would becone practically threshold insensitive. 
The turbulence intensity was never completely independent of the 
threshold setting since as the threshold increased the Macrodyne would 
only validate the larger particles which do not follow the flow as 
well. The threshold was set at the knee of the curve resulting in the 
optimum data rate, noise rejection, and particle sise validation. 
This corresponded to setting the threshold minimum validated sinewave 
amplitude to 90Z of the noise amplitude at the **Filtered Output'*. The 
noise level versus threshold setting calibration curve for this 
criterion is given in Figure 2.14. 

The overload discriminator setting did not make a significant 
improvonent in the turbulence intensity measurement. Large particle 
rejection was not important and the overload level was simply set well 
above the input signal levels. 

As mentioned earlier the velocity PDF was calculated for each 
data record. The PDF of the burst period was also displayed on the 
CRT terminal for each 2048 point record to check the Macrodyne 
opi^ration. The Macrodyne only accepted a window of Doppler burst 
frequencies determined by the input filters and counter range. This 
window had to be larger than the burst frequency bandwidth. Also the 
Bragg cell shift had to be high enough to obtain nonsero, unambiguous 
velocity measurements. All these criterion could easily be examined 
by viewing the burst period PDF. The boundaries of the PDF were 
displayed in both time and frequency units. The maximum and minimum 
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burst frequencies were checked egainst the input filters to prevent 
cutoff. The Bragg cell shift was adjusted so that the ainiaa Doppler 
frequency was above the ainiaun filter setting of 0.5 MBs. Display of 
the aasiaua and ainiatai burst periods allowed the operator to set the 
counter range to the highest sensitivity irithout overflowing the 
counter register. 

Finally the LDV velocity profiles were coapared with the hotwire 
aeasureo.nt: b . Also the velocity profiles were integrated and the 
integrated values %»ere checked against the aass flowrate 
aeasurenents. This trill be discussed in the next chapter. 

2.5.6 Particle Seeding . Cyclone aerosol generators described by 
Glass and Kennedy (1977) were used to suspend alumina (Al20j) 
particles into the flow (See Figure 2.16). Air forced thru a noxxle, 
mounted tangentially to a plexiglass cylindar, swirled the air and 
entrained the alumina particles at the bottom. The particle laden air 
was drawn out the to... Two generators %iere mounted in parallel and 
their exhaust lines were connected to a manifold. The manifold 
injected the seeded air downstream of the two sonic air noxzles 
(Figure 2.17b). To minimize contamination of the quartz windows, the 
alumina particles were injected into the air flow only during the 
acquisition of the LDV velocity-time record. A three-way solenoid 
valve directed the seeder secondary air either through the seeders or 
a bypass so that the total air flowrate was constant at all times 
(Figure 2.17a). 

The almnina particles (0.05 micron dia., 3.7 x 10 Rg/M , 
Linde Div., Union Carbide) were dried at 500 K for eight hours before 
each experiment to limit particle agglomeration. The entraining air 
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was also dried to prevent additional water abaorpti<m. The S.ld 
eicron fringe spacing prohibited particles over 1.5 wicriMis froa 
generating a Doppler frequency. Heasurewents by Glasa and Kennedy 
(1977) indicate that particle aggloweration occurs increasing the 
average particle diaweter by about a factor of 3. This would indicate 
an average particle diaweter of abou*' 0.2 aicrons was produced . One 
could safely assiae that the particle diaweter was below 1 wicron. 

2.5.7 Sources of Error . Errors in the LDV weasurewent can reault 
from any of the different facets of the weasurewent systew: optical 
alignwent, signal processing, gas wedita, and seed particles. 

Processor errors . A common source of error in frequency 
measurement is Doppler ambiguity broadening or transit time broadening 
(George and Lumley, 1973). The precision uhich the frequency of a 
waveform of finite length can be measured is limited. For a processor 
working in frequency space, a wave form N cyles long is effectively 
Fourier transformed and the relative width is of order 1/N. This type 
of broadening was avoided by using a period counter %ihich works in the 
time domain (Self and Whitelaw, 1976). The waveform of an individual 
burst can be '.imed to a precision only limited by the clock 
resolution. As discussed earlier (Section 2.5.3) this results in a 
0.2Z uncertainty. 

Optical errors . Besides the frequency measurement, the error in 
the velocity is dependent on the uncertainty in the beam crossing 
angle, 6 , the laser wavelength, A , and the Bragg cell shift .Vo. 
(Equations 2.1 and 2.6). The single sample uncertainty for the 
velocity oteasurement was estimated to be 0.4t. 

Fringe gradient broadening . Additional broadening in the Doppler 
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sigc :1 can result fron curved fringes in the focal voliae. Laser 
light propagates as Gaussian beans and the nintHnni vaista of the beans 
aiay not be at same position as the intersection point of the beams 
(Durst and Stevenson, 1977, Hanson, 1973). If the waists of the beass 
are not at the intersection poiut, the fringes %ra.ll be curved, the 
fringe spacing, d^ , will not be constant, and broadening of the 
Doppler signal will occur. The broadening was estimated to be O.IZ 
and will add to the ras turbulence. This apparent turbulence will 
only be important in the free stream above the layer. The gradient in 
the fringe spacing is too small to produce any discernible error in 
the mean velocity measurement. 

Velocity gradient broadening . Gradient broadening due to velocity 
gradients in the finite sized focal voluse can result in errors in the 
calculation of the mean and ms velocities. For individual 
realizations. Kreid (1974) approximated the error in the mean velocity 
to. 


u 





+ • • • 


2.11 


where u is the measured mean, ua is the actual aiean, 
probe volume diameter in the gradient direction, y. 
and u* ' are evaluated with respect to y in the center of the probe 
volume . 

For the streamwise velocity, at any given distance detmatream of 
the step, the maximum gradient will occur in the center of the shear 
layer. The velocity variation is linear and only the second term of 


and by is the 
The gradients u* 


equation 2.11 will be Important. The most severe gradients are near 
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Che edge of Che eCep end ere ■exisieed eC Che higheeC inleC velocicy 
(Oo ■ 22 M/e>. 

For exemple, ehen Do ■ 22 N/e eC 0.5 eCep heighCe dowieCreeB, Che 
■exiaum gredienc in Che leyer wee 1.3 x 10^ /e. The probe volute hec e 
dieaeCer of 96 aicrone in Che direcCion of Che gredienc. A aeen 
velociCy of 10 M/e will be 0.4Z in error. ISie gredience downcCre« of 
chie region were auch leee. AC one eCep heighc downeCreea, e gredienc 
of 4.5 X 10^/e wee Cypicel which givee a 0.05Z error eC e local 
velocicy of 10 M/e. 

CurvaCure effecce ere iaporcenc ec Che edge of Che ehear layer. 
The etrongeet curvature will be at cbe edge of the boundary layer at 
eeparation. The maximum effect of the firat tera of Bqn. 2.11 can be 
eeciaated there. The boundary layer profile ia laainar (See Section 
3.1). From a Blasius fit of the boundary layer at the highest Reynolds 
nunber (Rej^ * 3.7 x lO**, Uo • 22.2 M/s, u^ • 20 M/s) the aaxiaun 
curvature is -4.6 x lO’ /(Ms) «diere u' ■ 1.8 x 10** /s. This will 
produce a O.lZ error in the mean velocity for the worst case. 

The rus turbulence will also be affected by gradient broadening. 
Following Durst et el (1976), one assunes that the variation of the 
velocity in the probe volune can be expanded into a Taylor series. 
Since the laser beam has a Gaussian intensity profile, the velocity 
variation is weighted by a Gaussian probability function. The rms 
turbulence is Chen, 


h 

(u - Oa)* 



^ 6^ N** 



2.12 


The curvature term is small throughout the flow field. However, the 
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first tens in Bqn. 2.12 is iaportsnt in the boundary layer at 
separation and in the shear layer. The gradient broadening was 
subtracted in the boundary layer by calculating u* from the Tl-velocity 
data. A polynomial fit of three neighboring velocity points was used 
to determine the gradient. The maximum gradient broadening in the 
boundary layer ranged from Uj^/Uo “ 0.04 to 0.07 for the three 
Reynolds numbers studied. This amounted to about half of the 
turbulence measured in the boundary layer. 

Calculating u* directly from the data added noise to the 

turbulence intensity measurements. Therefore in the shear layer, the 
iT-*velocity data was fit with an error function by a least squares fit 
routine. The gradient was calculated from the error function fit and 
the broadening was removed. The maximum correction in the shear layer 

was about u^_/Uo * 0.03, 0.02, 0.005 for x/H “ 0.5, 1, 2, 

nns 

respectively. 

The skewness and flatness factors %Aiich are the third and fourth 
moments of the velocity PDF were also determined in the shear layer. 
Gradient broadening will also effect their evaluation. Again 
expanding the velocity in the probe volume into a Taylor series and 
averaging over a Gaussian probability one finds the skewness 
broadening is given by, 

(u-Ua)’ - ^by"u’*u" 


snd the flstness brosdening is given by, 



33 


(u-u-)' 


.Ik ^ 

16 ^“ + j_2g 


by*u’*u’'* 


4096 **y “ 


2.14 


Both terns in the skewness broadening were found to be negligable. 
The first term in the flatness broadening is not negligible. However 
if one normalizes the flatness factor by the uncorrected rms 
turbulence, the broadening due to the first term is removed. The rest 
of the terms in the flatness broadening are insignificant. 

Velocity biasing . In turbulent flows, Mclaughlip and Tiedeman 
(1973) first recognized that the mean and rms velocity averages will 
be biased because at higher velocities sK>re fluid and hence more 
particles are convected through the focal volume than at lower 
velocities. Time ave ge computations proposed by Dimotakis (1976) to 
circumvent the McLaughlin-Tiederman bias in high data rate regions has 
been presented in Section 2.5.4. The time weighted averages are 
correct for regions where the data rate is high. In the low data rate 
regions where time averages are no longer valid, simultaneous 
measurements of at least two and sometimes three components of the 
velocity must be made to properly correct for the biasing (Dimotakis, 
1976). A ringle component LDV system was employed, prohibiting this 
type of correction. In the low data rate regions, the nunerical 
averages (Equations 2.7 and 2.8) were reported without corrections. 

Bragg cell biasing . Velocity biasing is also caused by frequency 
shifting. The fringes move at a constant velocity and the number of 
realizations per particle will depend on the velocity of the particle 
(Meyers and Clemnons, 1979). In the limit, a stationary particle in 
the focal volume will produce an infinite number of realizations. 
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This velocity biesing effect will serve to cancel the 
McLsughlin-TiedenBan bias. The inportance of velocity biasing in the 
low data rate regions can be partially checked by exaaing the syosaetry 
of the velocity PDFs which were stored at each measurement location. 

Particle Considerations . The ability of the particles to follow 
the flow is an important consideration. The upper limit on the 
particle sise is determined by the particle inertia and the lower 
limit by Brownian motion, molecular slip, and the particle scattering 
cross section. The particle concentration must not be too high to 
affect the flow. Also the size distribution must be narrow since the 
velocity lag is dependent on the particle size. 

The cyclone particle generators using 0.05 micron alumina 
particles produced a particle diameter of about 0.2 microns (See 
Section 2.5.6). Mazumder and Kin^ch (1975) set the lower limit of 
particle diameter determined by siolecular slip at 0.1 microns. 
Brownian motion need only be considered for very low velocity (^0.05 
M/s) laminar flows (Durst et el, 1976). The particles were large 
enough to produce data rates up to 10 kHz in the free stream. Thus 
the particle size was not too small. 

The particle concentration was about 1 x 10*® /M* in the free 
stream. This represents an average separation distance of 2000 
diameters for 0.2 micron particles, causing a negligiole effect on the 
flow. 

The upper limit on the particle size is determined by the 
required frequency response The particles need to follow flow 
oscillations in the LaGrangian (moving) frame rather than the Eulerian 
(laboratory) frame which relaxes the size limitations considerably. 
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Using the particle flow solution by Hjelafelt and Mockros (1966) » 
Durst et el (1976) calculated the frequency response for various types 
of particles in both air flows and flaaes. Iheir calculations 
indicate that alumina particles 0.4 microns in diameter will follow 
flow oscillations (in the Bulerian frame) up to 10 kHs with IZ 
velocity lag in room temperature air flows. Particles of 0.8 micron 
size will provide the same performance in flmses. Thus 0.2 micron 
particles used here will have a frequency response in excess of 10 
kHz. 

The size distribution of the alumina particles was not measured. 
The boundaries of the particle size PDF were determined by the initial 
particle size (0.05 micron) and the fringe spacing (3.1 micron). Thus 
realizations could only occur from particles ranging from 0.05 to 1.5 
microns. Nazumder and Rirsch (1975) noted that the velocity lag 
varies with the size distribution and a %n.de distribution of particles 
might broaden the Doppler frequency. This may occur for highly 
accelerating flows with low turbulence levels. The flow field behind 
the rearward-facing step had turbulence levels exceeding 4Z everywhere 
and turbulence levels greater than 50Z in the shear laye^*. The large 
particle velocity lag was at most and would tend to narrow the 
Doppler spectrum due to the low frequency response of the larger 
particles. This is clear from the effect of threshold setting on the 
turbulence level shown in Figure 2.15. As the threshold setting is 
increased, only the larger particles are validated, and the turbulence 
level decreases. 

Combustion effects . Use of LDV in a combusting meditzs generates 
new sources of imprecision. Temperature and concentration variations 
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lead to index of refraction variationa that cauae movement and 
diaperaion of the laaer light. The probe voluae position fluctuates 
and the beams may be deflected such that they no longer cross. 
Propagation of the Gaussian beams through the index of refraction 
ihhomogeneities increases the waist diameter, distorts the fringes, 
and broadens the Doppler signal. These problems and others associated 
wit> combustion systems have been discussed by Self and Whitelaw 
(1976) and to some extent by Durst et el (1976). 

Beam deflections and the ensuing beam uncrossings did not appear 
to be a problem in the t%K>-dimensional combustor. The data rates for 
isothermal and reactive flows %iere essentially the same. The light 
deflections in the schlieren system (Figure 2.6) indicate that beam 
deflections of the order of 0.2 mm could have taken place in the 
burning shear layei'. Deflections of this magnitude would certainly 
result in some signal dropout. However, the probability of such beam 
uncrossings lujst be small because no significant signal loss was 
recorded . 

Durst and Kleine (1973) reported velocity biasing due to particle 
density variations in combusting flows. The high temperature, low 
density regions contain lower particle densities than the unburned 
regions. Uncorrected averages will be biased toward the velocity in 
the cooler turbulent eddies. At high data rates this bias was 
corrected by storage of che velocity-time record and time averaging 
(Equation 2.9 and 2.10). At low data rates, this biasing added 
additional uncertainty to the numerical averages (Equations 2.7 and 
2 . 8 ). 


Durst and Kleine (1973) also studied the effect of metal oxide 
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particles on the flame ve'^ocity. They found no significant change in 
the flame velocity in a premixed, stoichiometric, natural gas flame 
when using HgO particles at a concentration of 2 x 10^ /m’ . Thus the 
alunina particles used in this experiment at a maximua concentration 
of 1 X 10“ /M* had no influence on the combustion process. 
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Table 2.1 Conbaster Data Record 


Input 

Computer 

Teminal 


(TA.PA) Anient air temperature and preaaure 
(TT,PT) Teat aection temperature and preaaure 
(TF,PF) Fuel rotameter teaiperature and preaaure 
tranaducer aetting 
(TS,PS,T8) LDV aeeder rotometer temperature, preaaure, 
and diaplacement 

(TN,DN) Main air noszle temperaturea and diametera 


Deaired 

Conditiona 

Output 

Settings 


f(Uo) 

[(Phi) 

f(PHl,PN2) 

L(YF) 


Reference velocity 
Equivalence ratio 

Main air nozzle transducer settings 
Fuel rotometer displacement 


Output (MA) 

Conditions (MAT) 

.(MF) 


Main air mass flow rate 
Total air mass flow rate 
Fuel mass flow rate 


Table 2.2 Schlieren Optics Specifications 


Camera 

Test section 
view 

*1 

(mm) 


01 

(deg) 

02 

(deg) 

Movie 

Full 

220 

120 

3.2 

3.8 

Movie 

Enlarged 

283 

234 

3.2 

3.8 

Movie 

Super enlarged 

163 

363 

3.2 

3.8 

Still 

Full 

120 

1000 

3.2 

3.7 
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Table 2.3 LDV Data Record 


Conputer 
Acceae : 

Parallel FBurat counter range (4-7) 

Line icoaparator aetting (5/8) 

Interface [comparator accuracy (2) 

Record length (2048) 

Computer Fringe spacing (3.114 microns) 

Terminal Bragg cell shift (0.01 - 10 Mbs) 

Probe location (X,T,Z) 

Parallel plurst period Data Period 

Line 

Interface 

(x 2048 points) 
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Figure 2.5 Schlleren optical layout (dimensions In 
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Figure 2.7 Hot wire calibration curve. 
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Figure 2.8 Typical LDV dual beam optics. 
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Figure 2.9 Two-dimensional combustor LDV optics. 
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Figure 2.12 Macrodyne counter front end detector. 
(Photo cou-tesy of Macrodyne, Inc., Rotterdam, N.Y. 
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M^ure 2.13 r^acrodyrie counter lOQic and output 
(Photo courtesy of ?lacrodyne, Inc., Potterdam, 









Ini«n«tiy C%^ Dai 






OK QUALITY 



Figure 2.16 LDV cyclone particle seeders 
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Figure 2.17a LDV seeder air supply. 
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Figure 2.17b LDV seeder layout. 
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Chapter 3 

nSULTS AKD Discussion 


Three Reynolds number conditions (R^ * K3 x 10 » 2.2 x 10 » 
3.7 X 10*^) are studied for both reacting and non-reacting flow. The 
step height is fixed at 25 mm. In the reacting flow premixed 
propane-air is burned at an equivalence ratio of 0.57. The Reynolds 
number limits are determined by the stable burning modes of the 
combustor. All the inlet conditions are at room temperature. 

Complete profiles of the mean and rms averages of the u-velocity 
component are recorded for each flow condition. The velocity PDF is 
stored at most locations. The reacting flows are visual ixed with high 
speed schlieren photography. Measurements of the v-velocity component 
are made for one Reynolds number (R^ * 2.2 x 10**). Finally, the 
frequency spectra of the v-component velocity are obtained to assess 
the local large scale structure passing frequency at two inlet 
Reynolds numbers » 1.5 x 10**, 2.2 x 10**). 

3 . 1 Initial Coadit ions 

The state of the entry flow can have substantial effects on the 
behavior of mixing layers and reattaching flows. The condition of the 
boundary layer at separation can affect the shear layer's virtual 
origin, initial breakdown, groi^h rate, reattachment point, and 
turbulence development (see Bradshaw, 1966, Browand and Latigo, 1978, 
Batt, 1973, and Eaton, 1980a). Other variables such as test section 
geometry (Dimotakis and Brown, 1976), and free stream turbulence 
IChandrsuda et el, 1978) can also be important. In rearward-facing 
step flows, the major scaling parameters assessing rhe state of the 
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boundary layer are the aomentum thickness Reynolds number, Uq6/v, the 

boundary layer thickness, 5/B, and possibly the shape factor, 
* 

»12 /0. Here 6 1 b the 991 boundary layer thicknesB. The 

displacement and momentum thicknesses are £* and 0, respectively. 

Birch (1977, 1980) suggests that the plane mixing layer flow 
becomes independent of the initial conditions at very high Reynolds 
number (Re > 2 x 10^). Most investigations in mixing layers and 
backward- facing step flows are well below this value and initial 
conditions are important. The effect of the state of the boundary 
layer (laminar or turbulent) on the shear layer growth and development 
for Re < / X 10^ is still not wll understood. Batt (1975) and 
Wygnanski and Fielder (1970) report an increase in the gro%rth rate 
when the boundary layer is tripped and Browand and Latigo (1978) 
report a decrease. In this study the maximum '^eynold? number, * 

before curvature effects become large is 2 x 10^ making the 
documentation of the initial conditions important. 

The entry velocity profiles on the centerline for the three 
Reynolds numbers studied are shown in Figures 3.1, 3.2, and 3.3. These 
are numerical averages of the velocity. The time integrated averages 
are identical except near the step wall %#here the data rates are too 
low for the time averages to be meaningful. The inlet velocity is flat 
to ♦ 5X of the average freestream velocity. The variations in the 
mean velocity are due to the small contraction ratio (2:1) and 
nonuniformities in the steel %rool packing. Integration of the 
centerline profiles yields mass flowrates that are 2X - 32 lover chan 
the flowmeter measurement s , The free stream turbulence levels are 


about 21 for isothermal flow and 32 for reacting flow. The noise level 
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for the turbulence intensity measurements in the reacting flow is 
higher. 

Hot wire anemometry measurements of the entry velocity profile 
are made to check the LDV operation. Velocity profiles for 
non-reacting flow measured by both techniques are shown in Figure 
3.4. The agreement betmen the measurements for the three entrance 
Reynolds numbers varies from 0.3Z to 2Z. This is %#ithin the accuracy 
of the measurements which is 1.5Z for hot wire and 0.4Z for the LDV 
(see Chapter 2). 

The boundary layer for the low Reynolds number case is plotted 
against the Blausis profile for 0*0 and 0.57 in Figure 3.5. The mean 
velocities follow the Blasius profile and the shape factor of the 
profiles is the Blasius value (H^^2 * 2.6). The turbulence levels are 
corrected for gradient broadening (see Section 2.5.7) which is 
significant in the boundary layer. The correct ion is proportional to 
the velocity gradient (Eqn. 2.12) which is not a smooth varying 
function and the noise level of u^ms increased. The turbulence 
level is high in the boundary layer (^4Z) but flat. The heat transfer 
back to the step in the reacting flows did not alter the boundary 
layer significantly. The medium and high Reynolds number flows are 
shown in Figures 3.6 and 3.7. The medium and high Reynolds number 
velocity profiles depart from the Blasius values but the shape factors 
of 2.5 and 2.4 are still much larger than the turbulent value of 1.3. 
The turbulence levels increase in the ’ >undary layer for the higher 
Reynolds number flows. 

The boundary layer parameters are summarired in Table 3.1 for all 
the run conditions. The influence of combustion on the boundary layer 
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state is minimal. Under all conditions the turbulence levels in the 
boundary layer are high (4 - 81) yet the velocity profiles are 
basically laminar. This boundary layer state is similar to the 
tripped layer reported by Batt (1975). He used a wire trip producing 
high turbulence levels with a Blasius profile. Here imperfections in 
the profiled step (probably at the leading edge) might have tripped 
the boundary layer. The momentum thickness Reynolds numbers are 
around the critical value of 160 but are well below the transitional 
value of 390. The boundary layer never becomes fully turbulent. The 
condition of the boundary layer for all three Reynolds numbers is best 
described as "transitional". 

The initial boundary layer state can also have important effects 
on the stabilization of the reacting shear layer. Figure 3.8 shows an 
enlargement of the mixing layers for the reacting flows. The dark 
object protruding from the top is the ignitor. The schlieren 
sensitivity is minimized in Figure 3.8b so that only the flame is 
shown. The flame stabilizes in the boundary layer just before 
separation. With too large a boundary layer the flame could propagate 
upstream eausinr flashback. At low flow rates (Re„ < 1.5 x 10**) the 
ignition point cf the flame is found to wander in the boundary layer 
causing a low frequency flapping of the reacting layer. As the 
boundary layer becomes thinner , the ignition point of the flame must 
move do%mstream, and the flare could destabilize and blow off. 

Hot wire spectra of the rms streamwise velocity component in the 
entry region were taken for the non^reacting flow. In the study by 
Ganji and Sawyer (1980), the premix section had no packing or 
screens. The hot wire spectra revealed many resonant peaks from 
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acoustic waves in the tunnel (see Figure 3.9). These ''tuned** flow 
oscillations could trigger instabilities in the shear layer. 
Installation of screens and packing in the premix section (Chapter 2) 
eliminated this resonant turbulence. The hot wire spectra in the 
entrance flow for this study is given in Figure 3.10. The spectra is 
smooth and void of acoustic resonances. As will be seen later » under 
combustion the LDV rms velocity spectra exhibit some resonant peaks 
due to acoustic waves in the test section. It will be shorn that they 
are not of sufficient magnitude to affect the layer. 

3.2 Shear Layer Development 

The rearward- facing step flow develops initially very much like a 
free shear layer with zero velocity on one side (half jet mixing 
layer). Before curvature effects become large ( 6^yR « 1 > where R is 
the radius of curvature of the dividing streamline) the growth rate 
and velocity profiles are similar to the half jet (Eaton and Johnston^ 
1980a). Unlike the half jet flows, the upper velocity, U^, does vary 
slightly with x and more importantly the velocity in the recirculation 
zone is not zero. Since the mixing layer growth rates are dependent 
on the velocity difference across the layer, one would expect an 
effect on the growth rate. Turbulence levels in the recirculation zone 
are also much higher than the half jet mixing layer. The influence of 
these higher turbulence levels in the recirculation zone on the 
turbulence development of the mixing layer are not %rell understood 
(Eaton and Johnston, 1980a). 

As the layer continues to grow, curvature effects become large 
and the analogy to a mixing layer no longer applies. The layer curves 
down toward the wall and eventually impinges on the wall at the 
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reattachment point. The shear layer behind a re arv^^rd* feeing step is 
often called a "reattaching" shear layer. The recirculation sone site 
depends to a great extent on the rate of growth of the reattaching 
shear layer into the recirculation zone. The growth of the upper 
boundary of the shear layer will determine the rate of entrainment of 
fresh fluid into the layer. 

As in mixing layers, the turbulent structure of the reattaching 
shear layer is determined by the formation and development of large 
scale structures. The initial laminar shear layer breaks down from a 
Relvin-Helmholtz instability (Sherman, 1976) and the large scale 
structures form. These large scale structures are primarily 
two-dimensional (Brown and Roshko, 1974) and grow by fluid entrainment 
and coalescence (Winant and Browand, 1974). The growth of the large 
scale structures mark the boundaries of the reattaching shear layer. 
Their growth effects the recirculation zone size and the rate of 
spread of the upper boundary of the shear layer into the free stream. 
The turbulence in the layer results from the passing of these 
two-dimensional structures which contain a smaller scale, 
three-dimensional framework ( Breidenthal , 1978, Konrad, 1976). 

In the reacting layer, the incoming flow is premixed reactants. 
The flame is stabilized at the edge of the step. Near the step, the 
reactants and hot products mix in a laminar layer and burn. Once the 
eddies form the burning takes place in the two-dimensional structures 
(Ganji and Sawyer, 1980). The eddies fold in reactants and hot 
products as they roll downstream. As we will see, the heat oxpans^on 
causes the shear layer defined by the mean velocity field to be 
shifted downward and the recirculation size shortened. The large 
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scale structures react and expand above the shear layer. The density 
above the shear layer is lowered as the upper flame boundary extends 
beyond the boundary of the shear layer. The coalescence rate of the 
structures is also lowered by heat release. These and other effects 
of reaction on the reattaching shear layer will be detailed in the 
next sections. 

3 • 3 Mean Velocity Flow Field 

3.3.1 Sreamwi^e Velocity Flow Field . A typical streamwise 
velocity profile is shown in Figure 3. II. Both the numerical and 
integrated mean values are plotted along with the corresponding data 
rate. Above the shear layer the data rates are fairly high and the two 
means are identical. In the shear layer there are some differences 
betmen the two means but the values match closely. Only in the 
recirculation zone, where the data rates are typically low, are 
significant differences apparent. Clearly the time integrated values 
are suspect here since the data rates are too low to truly follow the 
flow. The numerical averages are also subject to velocity biasing 
errors in this turbulent region. In all the mean profiles the only 
significant departure of the two mean values occurs in the 
recirculation zone. Since the numerical means are more appropriate 
for this region, and the two values are the same elsewhere, all the 
velocity profiles reported will be numerical averages. 

Integration of the centerline velocity profiles in the isothermal 
flow yield mass flow rate estimates that range from 8X to 23Z below 
the flowmeter measurements. Off-center profile measurements indicate 
that this flow deficit can be attributed to lower recirculation rates 
near the side walls. A centerline profile at x/H ■ 3 is compared to an 
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off-center profile taken halfway between the wall and centerline 
(z/H * 1.8) in Figure 3.12a. The average bulk velocity in the p. of ilea 
is 20Z low in the centerline and 31 low at z/H * 1.8. The lower 
recirculation rates near the wall increase the average bulk velocity. 
The velocity in the recirculation zone (y/H » -0.67) for the same x 
position is shown in Figure 3.12b confirming that the recirculation 
rates are lower avay from the centerline. A complete volume integral 
involving sequential scans across the tunnel is not made but the lower 
recirculation rate near the wall can easily account for the flow 
deficit on the centerline. 

The aspect ratio (tunnel width to step height) is 6.9. According 
to de Brederode and Bradshaw (1972) at aspect ratios greater than 10 
three-dimensional effects are small. Here three-dimensional effects 
are not negligible, but as seen from Figure 3,12 the flow is still 
predominantly two-dimensional. ..:hlieren pictures of the reacting 
shear layer (Figures 3.8 and 3.40) also demonstrate the 
two-dimensionality of the flow. 

The centerline strearowise velocity profiles of the reacting and 
non-reacting flow at the three entrance Reynolds numbers studied are 
shown in Figures 3.13 through 3.15. The dashed line shown in the 
reacting flow indicates the upper edge of the flame boundary measured 
from long exposure schlieren photographs such as Figure 3.44. The 
flame boundary is well above the shear layer given by the mean 
velocity profiles. 

For the non-reacting flow the position of the zero streamline if 
determined from the centerline profile by integrating from the top 
wall. The results indicate that curvature effects do not become 
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important until t/H “4 or • 0.6, vhere x^^ is the reattacl v.er « 
length. This is consistent with the results of Eaton ard Johnston 
(1980a). The three-dimensional effects are too strong to give accurate 
quantitive data of the zero stteamline position. 

The maximum reverse velocity at each x position in the 
recirculation zone for the non-reacting flow is given in Figure 
3.16a. The maximum reverse velocity is 0.33Uo which is higher than 
reported by other researchers. A review of backward- facing step flows 
by Eaton and Johnston (1980b) reports a typical reverse vtlocity of 
0.25UO. The higher value is probably due tj the low aspect ratio of 
this flow. As seen in Figure 3.12b the rev 3 rse velocity is 0.3Uo on 
the centerline and decreases to 0.2Uo at z/H “ 2.0. 

The maximum reverse velocities of the reacting flow are higher 
due to heat release. The maximum reverse velocities range from 0.39Uo 
to 0.48UO for the three Reynolds numbers. Since the flow is not 
adiabatic these depend on the amount of heat loss to the walls. 

The streamvise volumetric recirculation &c each x position is 
plotted in Figure 3.17. For the isothermal flow, the reversed flow 
per unit width is O.I 6 U 0 H. This compares with a value of O.O 8 U 0 H 
reported by Moss and Baker (1980). Moss and Baker's reverse velocity 
(0.2Uo) is lower than normal (0.25Uo) which accounts for some of the 
difference. Most of the difference is a result of the low aspect ratio 
used here which causes three-dimensionalities in the recirculation 
zone. For example, the centerline profile in Figure 3.12a has a 
recirculation rate of 0.17UoH and the off-center profile has a value 
of 0.09UoH. 

The volumetric recirculation of the reacting flows are given in 
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Figure 3.17b tor completeness* The mass recirculation rates are 
dependent on the local density* The volumetric recirculation values 
ere higher at 0*16UoH to '^*24UoH« The volumetric recirculation in the 
X direction drops off decisively to zero between x/H • 4 and 5 
illustrating the shorter recirculation aone size for the reacting 
flew. 

From the mean velocity pre files information on the shear layer 
boundaries, growth rate, fluid entrainment, and reattachment length 
can be deduced. The effects of the heat of reaction on these 
reattaching shear layer properties will be considered in the follo%ring 
sections . 

3.3.2 Growth Rate * To quantify the gro%rth of the layer one must 
consider its boundaries. One measure of the thickness of the layer is 
the vorticity thicRn'^ss, 6 » given by. 
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where 


AU = 


3.2 


and and U 2 ate the velocities in 
respectively, and is the aa * 

jet mixing layers (U^ * 0) AU 
rearward--facing step configuration 
velocity in the recirculation zone 


the upper and low^r streams , 
m slope in the layer. In half 
is identical to U^. In the 
varies due to the non-zero 
(See Figure 3.16a). Here 6^ is 


estimated using the maximum velocity in th recirculation zone near 
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the shear layer edge. 

The vorticity thickneaa for the isothermal flows are plotted in 
Figure 3.18. The two lower Reynolds number flows grow more slowly 
initially. The lower Reynolds number shear layers are initially more 
stable and the layer breaks down farther downstream. This is confirmed 
by eddy formation position measurements made previously from schlieren 
observations (Ganji and Sailer, 1960). All the growth curves have a 
fairly linear region and then drop off when wall interactions become 
important. Only the high Reynolds number data exhibit a truly linear 
growth and is fit with a slope of 6^ * 0.28 in Figure 3.18. Estimates 
of the gro%fth rate in the linear region vary from 0.27 to 0.2^. These 
growth rates are much higher than those reported previously (Batt» 
1977). Only the growth rates of tripped half jet mixing layers (Batt» 
1975, Wygnanski and Fiedler, 1970) approach this value (6^* • 0.23). 
The excess gro%rth can b** traced to the non-xero velocities in the 
recirculation tone. Mixing layers at high Reynolds numbers are self 
similar and their growth rite is dependent on the parameters, 

^ . 3-3 

s « / Pj 

X - ( Uj - U2 ) / ( Uj +Uj ) 

For isothermal flow vs ■ 1) the growth rate is linearly dependent on X 
according to (Brown and Roshko, 1974); 
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6 * • 0.18 X 3 5 

For half jet mixing layers ( X* 1) the gro%vth .ate ia 0.18. In 

reamard- facing step flo%ra the negative velocity in the recirculation 
xone causes the value of X to vary between 1.2 and 2.0 in the mixing 
layer region of the flow. An average value of 1.6 could easily 
account for high growth rate reported (5^* « 0.28). 

Half jet flows with tripped boundary le/ers grow slower at first 
and then overshoot the asymptotic growth rate (Husain and Hussain, 
1979, Birch, 1980). Overall the growth rate is high in the developing 
region and causes the virtual origin to appear downstream of the 

separation point. The isothermal mixing layers seem to follow this 
pattern here. All the virtual origins are do%fnstream of the step 
which suggests that the boundary layer is indeed tripped. The virtual 
origin for flow developing from a laminar boundary layer is upstream 
of the separation point (Husain and Hussain, 1979, Birch, 1980). 
According to Birch (1980) at high enough Reynolds numbers 

(Re^ > 2 X 10^) the growth rate reaches an asymptotic value and the 

virtual origin is downstream of the separation point regardless of the 
initial conditions. The mixing layer region of the flow in this %iork 
is clearly not well developed. The Reynolds number is too low 
- 2 X 10^), the growth rates are high, gnd the virtual origin 
is downstrei of the step. 

The vorticity thickness plots of the reacting layers are shown in 
Figure 3.19. The initial growth of the layer is delayed significantly 
over the isothermal case. Again this is a local Reynolds number 
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effect. The Reynolds ntaber varies like T . ftoa thensocouple 

■easurements in the layer (Ganji and Sawyer, 1980), there is about a 

factor of 4 increase in temperature, sfeich results in an eightfold 

decrease in Reynolds nmber. Thus the shear layer breakdown is 

delayed over the nonreacting case, reducing the growth rate. 

At x/B ■ 2, the vorticity thickness for Reg ■ 1.5 x l(f is much 

greater. This might be due, in part, to a low frequency instability 

of the reacting layer mentioned earlier. The ignition point which is 

located just upstream of the step expansion (see Figure 3.8b) is 

unstable for Reg « 1.3 x lO**. The ignition point oscillates in the 

streamwise direction causing the «dtole reacting la3fer to flap at a low 

frequency. At Re„ ■ 1.3 x lO**, the layer seems stable but under close 
n 

observation some movement of the ignition point is detected* This 
might result in a very slight low frequency flapping of the layer 
which gives the large value of 6^ at x/H * 2« 

Estimates of the gro%ith rate of the reacting shear layers range 
from 6^' ■ 0*26 to 0.29- Again only at Rejj ■ 3.7 x 10^ is the growth 
truly linear and the data is fit with 6 * • 0.29. A direct c ':rison 

of the vorticity thickness plots for the high Reynolds ntasber case is 
shown in Figure 3-20. Initially the growth of the reacting layer is 
suppressed due to the lomr Reynolds number in the reacting layer. 
Surprisingly, the groisth rates of the two layers are nearly identical 
at 0.2S for the non^reacting layer and 0.29 for the reacting layer. 

Brown and Roshko (1974) reported a 2St reduction in the growth 
rate for a half jet mixing layer (X * 1) at high density differences 
(p^/Pl « 1/7). No reduction in gro%rth rat-‘ vas found in this study. 
This suggests that the effect of the density difference across the 
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layer ia counter balanced by the effect of heat release and volisetric 
expansion. 

For a reacting shear layer, the vorticity thickness is only one 
sieasure of the layer boundary. The eorticity thickness defines the 
region of strong momentum transfer but the regions of conbustion and 
heat transfer can be quite different. As seen in Figures 3.13 - 3.15 
the flame boundary is veil above the edge of the shear layer. 
Schlieren pictures such as in Figure 3.40 show reaction taking place 
in the large scale structures propagating inside and above the shear 
layer. Thus combustion is not confined to the shear layer defined by 
the mean velocity flow field. 

3.3.3 Entrainment . From Figure 3.13 - 3.15, it is evident that 
the reacting shear layer undergoes a shift downward into the 
recirculation xone. The upper boundary of the shear layer is plotted 
in Figure 3.21a for Re^^ 3.7 x 10^. The growth of the reattaching 
shear layer into the upper stream is dramatically reduced. All the 
reacting layers showed this behavior. 

In the reacting flow, the flame extends above the upper boundary 
of the shear layer and the density above the shear layer is decreased. 
The upper streai!! velocity Uj^(x) is nearly identical for both the 
reacting and non-reacting flow. Thus the mass rate of entrainment of 
fluid from the upper stream into the shear lay<*r is reduced for the 
reacting flow 

3.3.4 Reattachment . The reattachment length is defined as the 
point in the separation region where there is no net flow reversal. 
The reattachment length is determined to a large extent by the growth 
of the mixing layer into the recirculation tone. The lo%rer boundary 
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of the layer is difficult to assess due to the very gradual decrease 
in the slope « The center of the shear layer can be defined by the 
position uhere the velocity equals AD/2. The centerline position of 
the reacting and non-reacting layers at Re^ * 3.7 x 10^ is presented 
in Figure 3.21b. The center of the reacting layer is shifted do%m«ard 
propagating toward the wall nore rapidly. 

The re attachment lengths can be estimated from the velocity flow 
field plots given in Figures 3.13 - 3.15. They are stsmsarixed in 
Table 3.2. Velocity profiles are taken in increments of H near the 
reattachment point making the values accurate to ^ 0.5 H. The 
re attachment lengths for the reacting layers are reduced by 20 - 30Z. 
This affirms the general shift of the layer toward the lower wall. 

The momentum thickness Reynolds nisaber of the boundary layer at 
separation, Re 0 > is also given in Table 3.2. The reattachment length 
varies with the boundary layer state and therefore Re 0 . Eaton and 
Johnston (1980a) report a strong increase of with Re 0 for a laminar 
boundary layer. The reattachment length peaks for a transitional 
boundary layer and decreases slightly to an asymptotic level for a 
tu'^bulent boundary layer state. The reattachment lengths do not show a 
strong dependence on Re 0 here indicating that the boundary layer state 
is transitional. 

The reattachment lengths for the reacting flow increase about 20Z 
as Re increases. For a laminar boundary layer state at separation, 
Eaton and Johnston (1980a) found a 30Z increase in the reattachment 
length with Re 0 . The high temperature in the initial shear layer 
region decreases the local Reynolds number and could have placed the 
functional variation of in the laminar region. 
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Hich conbustion, the reattacfaaient lengths are shorter and the 
flov relaxes back to channel flov sooner (Figures 3.13 * 3.13) « 
According to Bradshau and Wong (1972), in isotheraal flows complete 
relaxation occurs at a distance of 30 times the shear layer thickness 
at reattachment • This corresponds to 50 step heights downstream for 
this flow. 

3.3.5 Transverse Velocity Flow Field . The transverse velocity, v, 
is plotted for Re|| 2.2 x lO^* in Figure 3.22. The isothermal 
profiles vary with y as expected. Initially the dividing streamline 
is parallel to the y axis, and ^ is zero in the center of the layer* 
As the dividing streamline curves downward, v is no longer 
perpendicular to it and there is a net downward flow in the center of 
the layer. In the mixing layer v varies like y near the centerline as 
predicted by a simple gradient diffusion model (Gbrtler, 1942). 

The reacting layer profiles have a similar behavior except that 
heat release and density variation cause the profiles to be 
asymmetric. The center of the shear layer undergoes a shift downward 
consistent with the u results. 

3.4 Turbulence Development 

3.4.1 Turbulence Intensity . The streamwise turbulence intensity 
profiles for are shown in Figures 3.23 - 3.23. In both the reacting 
and non-reacting flows the regions of high turbulence are bounded by 
the shear layer, broadening as the shear layer grows downstream. Near 
reattachment, the turbulence levels decrease and the turbulence 
profiles begin to take on the characteristics of turbulent channel 
flow. In the reacting flow, the reattachment lengths are shorter and 
more of the post-reattachment region can been seen. The domstream 
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profiles are ragged indicating that the 2048 point records are 
insufficient to produce a true werage. 

Gradient broadening of the turbulence intensity is removed in the 

shear layer (See Section 2.5.7). The velocity gradient is calculated 

by fitting an error function to the mean velocity data in the shear 

layer with a least squares fit routine. The broadening in the top and 

bottom boundary layers is not removed. The maximum correction in the 

shear layer at du/dyl is about u /Uo * 0.03, 0.02, 0.005 for 

max rms 

x/H * 0.5, 1.0, 2.0, respectively. 

The exact position of the maximum value of the turbulent 
intensity is difficul*' to determine but the general trend can be seen. 
In the isothermal flow the position of the maximum initially coincides 
with the centerline (y/H * 0). It then dips slightly toward the 

recirculation zone and rises back toward the centerline near 
reatrachment . These same trends have been reported previously (Eaton 
and Johnston, 1980a). The reacting layer behaves the same way except 
the dip is much more dramatic. The maximum value follows the shifted 
shear approaching the bottom wall much more closely and then rises 
toward the centerline after reattachment. 

The maximum turbulence value, %^s/Vo||nax> plotted in Figure 
3.26a for the non-reacting flow. The initial increase of the 
turbulence intensity is due to the formation of the large scale 
structures %ihich form earlier for the high Reynolds number flows. The 
Re Iv in-Helmholtz instability that results in the formation of the 
large structures is dependent on a critical Reynolds number based on 
X. The instability occurs first for the high Reynolds number flows. 
These results are confirmed by eddy formation position data dedviccd 
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from schlieren observations made in non-reacting flow by Ganji and 
Sawyer (1980). They found that the eddies formed at z/H « 0.63, 0.S5, 
0.47 for Reg* 1.5 x 10 , 2.2 x 10 , and 3.7 x 10 , respectively. 

For all the non-reacting experiments, the turbulence intensity 
reaches a naximun 1 to 2 step heights before reattachnent. The 
turbulence level for the low and medium Reynolds maber cases 
gradually increase to this maximum. This monotomic increase of the 
turbulent intensity to a maximum 1 to 2 step heights before 
reattachment has been found in other rearward-facing step experiments 
such as Rim et al (1978), Eaton and Johnston (1980a), and Etheridge 
and Kemp (1978). A complete summary of backward- facing step turbulence 
levels has been given by Eaton and Johnston (1980b). 

The turbulence development of the high Reynolds number flow does 
not reach a single peak. The levels peak early at x/H •' 0.5 and then 
decay and reach a second maximum at x/H » 4. Bradshaw (1966) found a 
similar behavior in a mixing layer formed at the edge of a nozzle. 
Both the turbulent shear stress and u-turbulent intensity were 
double-peaked for the mixing layer formed from a laminar boundary 
layer. Bradshaw suggests that this second maximum "marks the 
establishment of the shear-producing part of the turbulence spectrum". 

4 

Examining the full turbulent profiles for Rej| » 3.7 x 10 in Figure 
3.25a» one notices that the profiles at x/H * 1-3 are flat topped and 
even double-peaked in the shear layer. The rest of the data is 
single-peaked which make these profiles suspect. Thus more 
measurements need to be made to confirm whether the turbulence growth 
is double-peaked and whether the second peak marks the emergence of 
fully three-dimensional flow. 
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nte naximuD turbulence level in the eheer layer for all three 
Reynolds aumbers is u^^g/Uo ■ 0.28. Turbulence intensity values of 
0.21 are sK>re typical for backward- facing atep flows (Baton and 
Johnston 1980a). A number of conditions lead to this high level. The 
freestream turbulence is high at 2-3X. Also the aspect ratio is low 
(6.9) which causes significant three-dimensional effects and results 
in a high reverse velocity. The velocity difference across the layer, 
AU, is 1.33UO as opposed to 1.2Uo for typical flows. Renormalising the 
turbulence levels would decrease the maximums by lOZ. 

Until curvature effects are strong, the shear layer formed at a 

rearward- facing step is very similar to a plane mixing layer. However, 

the maxiffliBD turbulence level found in backward-facing step flo%rs 

(urjQg * O.2IU0) is normally higher than plane mixing layers 

(uj^g ■ 0.171'o, Batt, 1977). The higher turbulence levels for 

backward-facing steps is a result of the reentrainment of pressure 

driven, highly turbulent fluid into the recirculation zone (Bradshaw 

and Wong, 1972). For half jet mixing layers the fluid on one side is 

at rest and the turbulent intensities are normalised by the constant 

velocity difference across the layer AU ■ Uo. In the backward- facing 

step flows AU is not constant and is typically 1.2Uo (AU • 1.3Uo in 

this experiment). Renormalising by the local value of AU «rould 

decrease the turbulence levels. Not only is the fluid in the 

reciculation zone not at rest, but is highly turbulent (u ■ O.lUo) 

nns 

which also increases the turbulences levels in the shear layer. 

The turbulence levels for all three flows decrease rapidly 1-2 
step heights before reattachment . This has been reported by other 
researchers (Eaton and Johnston, 1980a). The tt/o explanations proposed 
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to explain this decrease are shear layer curvature and splitting of 
the eddies near reattachment (Eaton and Johnston, 1980a). Castro and 
Bradshaw (1976) found curvature in shear layers to be stabilizing and 
cause a rapid decay of turbulence. The curvature of the layer is 
strong in the reattachment region which decreases the turbulence 
levels. Bradshaw and Wong (1972) reported a splitting of the eddies 
near reattachment, calculating that the length scales decreased by a 
factor of two. These eddies produce most of the shear stress and their 
splitting results in a decrease of the turbulence near reattachment. 

In this study, measurements of the eddy frequency (Section 3.3) 
show no reduction of the eddy wavelength in the reattachment region. 
The eddy structures as seen from the LDV and schlieren measurements 
maintain their integrity throughout the reattachment zone. Thus the 
primary mode for the turbulence decay is a result of the curvature of 
the shear layer. 

The reacting flow (Figure 3.26b) has higher turbulent intensity 
values of 30-33Z as a result of volimietric expansion in the layer. The 
data exhibits two separate growth and decay periods suggesting at 
least four mechanisms controlling the development of the turbulence. 
The first growth period is delayed over the isothermal flow as the 
large scale structures form later (Table 3.3). The subsequent 
turbulence decay is connected with the two-dimensional structures 
burning, coalescing, and growing as they are convected through the 
probe volume. The second growth region probably results from the 
introduction of small scale three-dimensional turbulence, Konrad 
(1976) found a three-dimeusional instability in the large scale 
structures formed in a plane mixing layer that increased turbulent 
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mixing by 25Z. Since the turbulent Schmidt number should be near unity 
this would also increase the turbulence intensity. For the three 
Reynolds numbers, the transition would occur at u/H ■ 1.3 * 2.5 for 
room temperature air. Heat release would delay the transition and the 
high turbulence in the boundary layer would tend to advance it. If the 
second growth period is due to Konrad *s instability, the Reynold's 
number should affect its location. There is not enough data in Figure 
3.26b to determine whether there is a Reynolds number effect. Again 
the turbulence decays as the curvature effects become strong, 1-2 step 
heights before reattachment. 

The turbulence intensity profiles of the transverse velocity are 

4 

given in Figure 3.27 for Re^^ * 2.2 x 10 . The values of vrms urms 

are the same in the freestream but is much less nearly everywhere 

else. In the recirculation zone at x/H * 2 and 5, is as large as 

rms 

The values of v /Uol are given in Figure 3.26. The 

rms rms I max ^ ^ 

transverse turbulence does not grow significantly until well after the 
large structures are formed. The isothermal values are about 30Z below 
the streamwise turbulence data but the maximum occur at the same 

position. These results are consistent with other backward- facing step 

data (i.e. Etheridge and Kemp, 1978, Kim, Kline, and Johnston, 1978). 
For the reacting flow at x/H » 0.5, the .Idies have not formed yet and 
there is no clear maximum value of in the shear layer. The value 

of v^g/Uo at the same y position as u^g/Uol is given in Figure 
3.26b. The doimstream values of v^g are similar to the streamwise 
data. The ^ which is about one step height beyond 

reattachment looks more developed than the u^g data indicating that 
the reattachment affects y^ms stroi^gly. 
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3.4.2 Probability Density Functions . Typical PDFs of the 
streamwise velocity component for the isothermal layer are displayed 
in Figure 3.28. The probability density^ P» on the vertical scale is 
unitless and defined such that the integral of (P)d(U/Uo) is 1. The 
weighting factor used to correct the LDV counter biasing discussed in 
Section 2.5.4 is not exactly correct and the PDFs still have 
systematic noise. To correct this problem the PDFs of the Doppler 
burst period should be stored and the PDF of the velocity calculated 
with velocity intervals corresponding to unit increments of the 
digital burst period. All PDFs from counter processors %rill contain 
systematic noise independent of sample size unless PDFs are calculated 
in this fashion. 

The positions of these Pdfs in the shear layer are indicated in 
Figure J.29c. The skevmess and flatness factors determined from the 
PDFs are displayed in Figures 3.29a and 3.29b. The skemess which is 
the third moment of the PDF indicates symmetry. The flatness factor or 
fourth moment describes how fast the distributic i decreases in the 
wings. A Gaussian distribution has zero ske%mess and a flatness 
factor of 3. The PDFs above and below the shear layer are Gaussian 
(”a” and ”f”). Near the edges of the shear layer (PDfs ”b”, ”d”, and 
”e'*) non-Gaussian behavior is seen. On the low velocity edge PDF ”b” 
is positively skewed with a long tail on the high velocity side. PDFs 
"d" and *'e" on the high velocity edge of the shear layer have tails on 
the low velocity side and «re negatively skewed. The flatness factors 
are high at these positions as the distribution falls off more slowly 
than a Gaussian 03). In the center of the layer (PDF ”c”) the 
turbulence is the hignest but the distribution is Gaussian. All the 
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skewness and flatness factors are normalised by the turbulence 
intensity calculated from the PDFs. 

The PDFs shown here are not corrected for velocity gradient 
broadening which will increase the width of the PDFs near the step* 
Gradient broadening will also affect the skewness and flatness factors 
and its contribution to these factors is given by Equations 2*13 and 
2.14. Gradient broadening of the skewness is estimated from an error 
function fit of the u-velccity data. The maximum contribution 
(Re^ * 3.7 X 10**, x/H = 0.5) is uVu^^ ■ 0.04 which is negligible. 
Gradient broadenin;.; of the flatness factor is removed by normalizing 
the flatness factor by the gradient broadened turbulence intensity. 
All the values of this section are calculated from the 
PDFs and are gradient broadened near the step. The skewness is also 
normalized by this value of u^^g. In the center of the shear layer at 
x/H ^ 0.3, this causes a normalization error giving values of u^ wh* !i 
are smaller by factors of 1.8 to 4.1. Beyond one step height the 
normalization error is negligible. The normalization error near the 
step does not change the sign of i*^ only its absolute magnitude. The 
character of the variation of the skewness in the layer is unchanged. 

The non--Gau8sian behavior of the turbulence near the ed^es of the 
shear layer in Figures 3.28 and 3.29 can be explained by considering 
the large scale structures. They form quasi*periodically at 
x/H ^0.3. The large structures intermittently grow to a sufficient 
size to affect the edges of the shear layer. When this happens th<v 
convect low velocity fluid to the upper edge (**d**) and velocir.i* 
fluid to the lower edge (**b**). The velocity of ♦♦ '^luid convecled by 


the large scale structures departs more radically from the mean 
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velocity than the velocity variation due to aaall scale turbulence and 
the tails are produced. The variation of skewness and flatness shown 
in Figures 3.29a and 3.29b is consistent with data in isothensal free 
shear layers (Spencer and Jones, 1971) 

The skewness and flatness factors for the low Reynolds nuBiber 
flows are given in Figures 3.30 and 3*31. The non-^reacting flow is 
negatively skewed in the upper edge of the shear layer and positively 
ske**^d in the lower edge. This behavior maintains itself even after 
reattacbcent . The non-Gaussian behavior is much stronger in the upper 
part of the layer where the gr^^^ents are stronger and the turbulence 
is lower. The flatness in the high velocity edge is much 
greater than 3. The flatness factor in ^'he recirculation rone is 
nearly Gaussian except nei^ the step (x/H < 1). 

The skewness of the turbulence in the reacting flow looks quite 
different first. However there are similarities id?en one considers 
that the reacting layer is shifted domward well below the centerline 
of the test section. Examining the u-velocity profiles in Figure 
3.13b, the "‘DFs ait positively skewed at the lower edge of jhe shear 
layer an:^ negatively skewed near the upper edge as in the isothermal 
layer . 

The non-Gaussian behavior for non-reacting and reacting flow are 
very similar at x/H * 1. Fir comparison the reacting flow PDFs are 
shuwn in Figure 3.32 at s Jiilar locations as ;$iven for the 
non-reactin^ PDFs in Figure 3.28. The position of ..he PDFs in the 
shear layer, the ske%/nes8 and flatness factors, anci the turbulenct: 
intensity are given in Fig ares 3.?3a ch. u 3.33d. The variation of the 
skewness and flatness factors is the * as the non-'reacting flow at 
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*/H - 1. 

There are two a«jor differences in the nonr^aussian behavior in 
t^e reacting flow. In Figure 3.30b, for x/H >1, tne skewness has a 
positive peak at the upper boundary of the shear layer given by the 
n-velocity prot les. The negative peak in the skewness is slightly 
below this boundary. To illustrate this behavior, the PDFs and their 
moments are displayed in Figures 3.34 and 3.35. As before PDF *V is 
skewed positively at the lower edge of the shear layer. However, just 
below the upper edge the skewness is negative C'c") and at the upper 
boundary the skewness is positive (”e"). 

The existence of this second positive peak is a result of 
combustion in the layer. As seen from schlieren observations, 
combustion occurs %rithin the large scale structures as they convect 
downstre^tm. The upper boundary on the gro%fth of these large scale 
structures can be determined from long exposure schlieren photographs 
such as Figure 3.44. This boundary represents the farthest excursion 
of the eddies into the upper inviscid stream and is the edge of the 
flame. The boundary as seen in Figure 3.13 is above the shear layer 
given by the "u-velocity profiles. The expansion of the eddies due to 
reaction fills out the upper part of the shear layer and feeds the 
momentum transfer across the shear layer. At Che same time the 
semi-periodic passing of the expanding fluid betwen the top of the 
shear layer .ind the upper flame boundary produces a high velocity tail 
on the Pi'Fs (See Figure 3.34, PDF "c”). The positive peak occurs 
between the upper edge of the flame and the cop of Che shear layer 
(See Figures 3.30b, 3.34 and 3.35). Thus the second positive skemess 
peak results from the impelling effect of the large scale structures 
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expanding between the top of the shear layer and the upper boundary of 
the flame. 

The flatoesB factor in the reacting layer is high from the top of 
the shear layer to the uppe*. thermal boundary even though the skewness 
changes sign (Figure 3 31b). For example, PDF in Figure 3.34 has 
nearly tero skewness but tails on both irings of the PDF produce a high 
flatness factor. At this position the large scale structures convect 
both high speed fluid from reaction and low speed fluid from the 
recirculation xone. Similar behavior in the skewness and flatness 
factors for reacting flow can be seen at the medium an'* high Reynolds 
numbers in Figures 3.36 and 3.37 except that the positive skemess 
peaks above the shear layer for the high Reynolds number flows are not 
as evident in the downstream post ions. 

The second major difference is near the step at x/“ ■ 0.5. In 
all the reacting flows the skevmess has negative peaks at both the 
lower and upper edges. The skewness has a positive peak near the 
center of the layer. This region is more clearly examined in Figures 
3.38 and 3.39. Typical PDFs are given in Figure 3.38. At the lower 
edge of the layer (**b") the PDFs are negatively ske%ied. The PDFs are 
positively skewed just below the center of the layer ("d"), Gaussian 
in the center (between "d'* and 'e*'), and negatively skewed near the 
top (”6**). The positively skewed PDFs near •M” are all bimodal. All 
PDFs at this location in the reacting flows are bimodal. Everywhere 
else the PDFs are e A-peaked. 

Two explanations for the double«>peaked nature of the PDFs near 
the step can be given. The intermittent formation of the large scale 
structures due to the instability of the shear layer could cause a 
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biaodal PDF. The foxMtion of the etrocturee pcoducee repid bvniiiig, 
high local heat releaae, and npanaion. The averse velocity of this 
•vent would be higher than the velocity in the laninar l^er when 
there is no breakdown. 

Secondly, the biaodal behavior could be a result of low frequency 
flapping. As discussed earlier, the reacting l^fer flaps slightly at 
the low Reynolds nuwber condition. Slight flapping of the thin shear 
layer at this position could produce a double-peaked velocity 
distribution. Since the bimodal PDFs are found for all three Reynolds 
nunbers and fron the schlieren aovies there does not seen to be any 
flapping at the higher Reynolds nusbers, the first ezplanatiou is wore 
plausible. 

3.5 Large Scale Structures 

The presence of large scale structures in free shear l^rers and 
in reacting shear layers has been firmly established (Brown and 
Roshko, 1974, Ganji and Sawyer, I960). The large scale structures 
form early in the shear layer from a Kelvin-Helmholts instability 
(Sherman, 1976). The structures are t%io-dimensional in nature and 
grow by coalescence and entrainment as they are convected through the 
layer. For premixed combustion in the layer, the reaction occurs in 
the eddies as they entrain premixed reactants and hot products. Heat 
release expands the structures and contributes to the growth of the 
l.'yer. 

Since the large scale structures control to a great extent the 
turbulence structure, a major goal in this study is to assess the 
effect of combustion on their evolution. The large scale structures 
are analysed with high speed schlieren photography and spectral 
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analysis of the U>/ signal. 

3.5.1 Schlieren Photographic Observations . Schlieren photographs 
of the reacting flow at the three entrance Reynolds nmbers are shorn 
in Figure 3.40. The exposure time varys from 49 to 74 ys. The large 
scale structures are prominent throughout the flow field. The rod 
shaped object protruding from the top is the ignitor. The dark bands 
highlight the high tmnperature gradients due to combustion. The 
combustion is taking place in large scale structures forming near the 
step. The fluid above and below the layer is fairly isothermal. 

The initial shear layer can be seen more clearly in Figure 3.S. 
The schlieren is desensitised in 3.8b and shows the point of ignition 
at the edge of the step. The large scale scale structures form as the 
initial shear layer breaks down. The structures grow by fluid 
entrainment, heat expansion, and coalescence. 

The mechanism of coalescence is illustrated in Figure 3.41. The 
sequence of photograghs are spaced 0.3 ms apart. The coalescence 
process which is the combining of two neighboring eddies is also 
called amalgamation or pairing. The pairing process is found in 
isothermal free shear layers (Winant and Browand . 1974) and in 
reacting shear layers (Ganji and Sawyer. 1980). The upstream eddy 
overtakes the do%mstream eddy and the two combine into one eddy with 
twice the original sire. Brown and Roshko (1974; reported eddy pairing 
to be one of the major gro%rth mechanisms in free shear layers. 

At the high Reynolds number flow, the large scale structures do 
not appear to be as ordered as the low Reynolds number flow (Figures 
3.8c and 3.42). This is probably due to three-dimensionalities in the 
breakdown of the layer near the step. For all the reacting flows the 
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flaae appears to atabilisa at aore than one region on the step edge 
allowing the layer to break down at different times across the flow. 

The high Reynolds nimber flow seeas aore susceptible to this 
three~diaensional breakdown causing the blurred pictures in Figures 
3.8c and 3.42. 

The average shedding position is detetained froa the schlieren 
movies. RistograiBs of the shedding position (Figure 3.43) from movies 
at Reg * 1.5 z lO" and 2.2 x 10^ are Gaussian. The shedding position 
decreases with Reynolds nusber as the shear layer is less stable at 
higher Reynolds numbers. 

Long exposure photographs are taken to deduce the flaae 
boundary. Figure 3.44 shows a schlieren ^otograph of the five where 
the bright areas indicate high temperature gradients. The top edge is 
the flame boundary and is uppermost propagation of the large scale 
structures into the upper free stream. This boundary al-rO represents 
the upper boundary of the thermal layer since the primary mode of heat 
transfer is convection of burning gases in tht large scale atructures. 

The flame boundary of the three Reynolds number flows is shown in 
Figure 3.45. The vertical scale is greatly expanded to emphasize the 
differences in the spreading rates. The flame spread decreases for 
the higher Reynolds numbers since the stoichiometric ratio is constant 
(<t> * 0.57). The flame spread for the low Reynolds nusd>er flow is 
higher than one would expect from the Reynolds number decrease. 
Slight low frequency flapping is observed at this condition idiich 
results in the higher flame spread. No flapping is observed at the 
other conditions. 

The average passing frequency of the eddies is measured at^^ 
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different dovnstream post ions. Using a aotion analyser a vertical 
line is positioned at a specified location on the aovie frame. Each 
time a eddy passes the line, the frame nuisher s recorded. The 
passing frequency is calculated as the inverse of the time between 
eddies. Histograms of the passing frequency are shown in Figure 3.47a 
and 3.48a. The histograms are positively skewed near the step and 
more Gaussian downstream. This passing frequency data is be compared 
with the passing frequency determined from LDV spectra in the next 
section. 

3.5.2 LDV Spectra . The frequency spectra of the vertical velocity 
compon nc is recorded to measure the large scale passing frequency in 
the layer. The structures are SK>st easily detected by positioning the 
LDV probe volume in the center of the shear layer and measuring the 
vertical velocity component. A spectrum taken in *^he isothermal layer 
is shown in Figure 3.46. The broad peak in the spectrum is generated 
by the large scale structure. 

The LDV spcctum for the reacting flow at the same position 
(Figure 3.47b, x/H * 2) has a broad peak superimposed with narrow 
acoustic peaks. The evolution of the large scale structure in the 
reacting flow is more cle «rly revealed in Figure 3.49. A broad peak 
in the spectra is first seen at x/H * 1 and the peak decreases in 
frequency with domstream distance. The width of the broad peak is 
about the same as found in the non-reacting flow where no acoustic 
evident. Narrow peaks in the spectrum from acoustic 
1 X..: -n the combustor are more pronounced near the step and 

Af j-p*itude as x/H increases. . 


jouscic resonances do .lot seem to affect the 
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dev«lops«nc of the large acele atuctorea. The broad peak 
correaponding to the large acele paaaing frequency appeara at x/H ■ 1 
and anoothly decreaaea in frequency through the layer. The broad peak 
ia alwaya nuch wider than the narrow aeoratic peaks indicating that 
the large acala atructnrea are never cowpletely coupled with the 
acouatica. The acouatic aodea are fixed in frequency * '^ughout the 
layer. 

When coupling of acouatica and the large acale developxent doea 
occur, the effects are draxatic. Keller et el (1981) excited the 
reacting ahear layer in thia apparatus with a loud speaker. Using a 
driving frequency of 200 Be at a relatively saiall power level, they 
were able to couple the large scale passing frequency with the 
acoustic frequency. High speed aovies of the reacting ahear layer 
disclosed that the large scale structures were completely coherent. A 
large increase of the growth rate took place idien the eddie passing 
frequency matched the acoustic driving frequency. The structures had a 
constant wavelength beyond this point and no coalescence is observed. 
Coupling of acoustic perturbations and the large scale stucture 
passing frequency has also been reported in mixing layers (Fielder et 
el, 1977). None of these coupling effects are observed here. 

The eddy formation posi'.. t;, Xg/H, is eatimated from the LDV 
spectra. The first appearance of a broad peak in the spectra is 
approximately the average eddy she Ming position. The estimated 
values are given in Table 3.3 along #ith the more accurate schlieren 
determinations. Measurements of the LDV spectra are taken every 1/2 
step height near the step and the shedding positions match within that 
uncertainty. Again the shedding position increases with decreasing 
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Reynolds ounber. In the reacting flow the local Reynolds number is 
decreased by a factor of eight from the iourfold increase in 
tmsperature. The decrease in Reynolds number increases for the 
reacting flow. 

The passing frequency from both the schlieren movies and LDV 
spectra are compared in Figures 3.47 and 3.48. Near the step for the 
medium Reynolds number flow (Figure 3.47) the schlieren sieasurements 
give a passing frequency 50Z higher than the passing frequency 
measured from the LDV. At x/H ~ 7 the difference is still 3SZ (Figure 
3.48). The Strouhal number, Str^ « Hf/Uc, is compared for the two 
measurements in Figure 3.30. The Strouhal number data from the 
schlieren movies does not have a similar shape for the two entrance 
Reynolds numbers and is consistently higher than the LDV data. 
Three-dimensional breakdown of the reacting layer probably causes some 
of the discrepancy. As mentioned before, the flame seems to stabilize 
at more than one location across the step. A close examination of the 
schlieren movies indicate that the layer breakdo%m is not always 
two-dimensional. Since the LDV measurements of the passing frequency 
are made at a single point , they are not effected by the 
three-dimensionalities. The second problem with the schlieren 
histograms is that they are not true spectra of the passing frequency. 
There is no correlation of motion with wavelengths longer than the 
distance between eddies. Thus the estimated passing frequency is 
higher than one determined from a true frequency analysis. 

The Strouhal number, Strj^, for both non-reacting and reacting 
flow are plotted in Figure 3.31. The eddies in the reacting tayer 
form later but at the same frequency, Stro* 2.0, as the isothermal 
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flow. Ihis •uggesca that the acoustic peaks lAich are visible only in 
the reacting flow do not affect the eddy formation frequency. The 
isothermal data appears to vary like a mixing layer Strg '^l/x between 
s/H *1-4. Specifically the Strouhal nusber should be given by 
(Brown and Roshko, 1974), 


Strfi 


Hf OH 
Oc (x - Xq) 


3.6 


where Xg is the virtual origin and a is a constant. In order to 
determine the virtual orgin. Equation 3.6 is rearranged as. 




Hf 


— ( — ) - — (^ 


3.7 


The isothermal data is plotted in that form is Figure 3.52. The values 
of a and are determined from a least square fit of the data between 
x/H ■ 1 4. The growth rate of the shear layer (Figure 3.18) is only 
linear in :hat region. The virtual origin in Figure 3.52 is downstream 
of the step as in the gro%rtk rate data suggesting that the boundary 
layer is tripped (S*e Section 3.3.2). The slope of the least squares 
fit is nearly identical for the two different entrance conditions. 

The Strouhal nianber, Str^, in the isothersial layer is plotted 
with respect to the virtual origin in Figure 3.53 for x/fi * 1 ‘■4. The 
Strouhal masher clearly has a l/(x - Xg) dependence given by. 


Str„ - 1.28H 

(x- Xo) 


3.8 


On the average everytime the distance from the virtual origin is 
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increased by 2*6 the Strouhal number is halved. 

In the isothermal layer coalescence seems to continue beyond 
x/H • 4 (Figure 3.51). There is not a clear trend in Strn for 
x/H > 5, as the convection velocity of the eddies, Uc, is difficult to 
estimate near reattachment. Bradshaw and Wong (1972) reported that the 
eddies are split in t%ro near reattachment. The rapid decrease in 
turbulence near reattachment is attributed to the ^bifurcation of the 
shear layer at reattachment**. They found at reattachment the length 
scale of the energy producing eddies is halved. As seen in Figure 3.51 
eddy splitting does not take place. The data suggests that the 
wavelength actually increases slightly through the reattachment 
region. Thus the rapid decrease in shear stress near reattachment is 
not a result of eddie splitting. 

The Strouhal number of the reacting data (Figure 3.51) does not 
show the simple l/(x - Xq) dependence in the linear growth region 
(x/H * 1 - 3, Figure 3.19). The Strouhal number only decreases 
slightly between x/H * 1 - 3. In Figure 3.54 the Strouhal number is 
replorted with respect to the virtual origin of the layer estimated 
from the growth rate data (Figure 3.19). Only the data in the linear 
growth region (x/H * 1 • 3) ic used in the curve fit. The coalescence 
rate is much lower than the isothermal layer having a weak dependence 
on the streamwise coordinate: 


Strn 


1.79H 

(x-Xo)®*“" 


3.9 


Thus the effect of combustion and volumetric expansion is to greatly 
reduce the pairing process as a growth mechanism in the shear layer. 
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The shear layer grows priaarlly by entraiiment of fluid and ▼oliatetric 
expansion as the fluid reacts. 

After the fairly constant region in the reacting shear layer, 
Str|| decreases fron x/H ■ 3 to x/H “ 7 (Figure 3.51). The eddy 
wavelength increases in this region froa heat expansion. The eddies 
are no longer retarded by the recirculation sone as they continue to 
bum and accelerate. 

The Biean spacing, 1, between the eddies in mixing layers is just 
1 ■ Uc/f. For the isothermal layer the spacing will be (Equation 
3.8): 


1 « 0.78 (x-Xq) 3.10 

This is a larger spacing than found in two-stream mixing layers but 
similar to a single stream mixing layer formed at a round jet. Brown 
and Roshko (1974) estimated the toean spacing to be 1 ■■ 0.92x from 
Rolpin's (1964) data. 

The mean spacing in the reacting layer (x/H * 1 - 3) is nearly 
constant (Equation 3.9), 

I - 0.56 (x-Xo)®-“’ 3.11 

emphasizing the low rate of coalescence. 

The mean eddy spacing in isothermal layers is linearly dependent 
on , 1 ■ ^ Brown end Roshko (1974) compiled mixing layer data and 
estimated the constant, c, to be about 3. This constant can be 
written in terms of the parsing frequency as. 
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c - 1/6^ . 

The inverse Strouhal number based on vorticity chickness, Uc/(f6J, is 
shown in Figure 3.55 for the non-reacting and reacting shear layer. 
The isothermal data approaches a value of c ■ 2.5, %«hich means the 
eddie mean spacing is 2.5 times the vorticity thickness. The reacting 
shear layer shows no constant value since the coalescence rate is very 
small and the wavelength, 1, is nearly constant %rith distance. 
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T^le 3.1 Batrj Flow CoDditiooa 


Uo 

(M/s) 

♦ 

Rsh 

**e 

6* 

(m) 

«12 

Free- 

Strv 

f w w 

Bound. 

Layer 

9.1 

0.0 

1.5 

X 

lo" 

140 

0.63 

2.58 

0.015 

0.040 

9.1 

0.57 

1.5 

X 

lO** 

150 

0.62 

2.50 

0.030 

0.045 

13.3 

0.0 

2.2 

X 

lo' 

180 

0.50 

2.46 

0.020 

0.050 

13.3 

0.57 

2.2 

X 

10* 

190 

C.55 

2.48 

0.0^0 

0.055 

22.2 

0.0 

3.7 

X 

10* 

280 

0.47 

2.42 

0 020 

0.080 

22.2 

0.57 

3.7 

X 

10* 

250 

0.42 

2.41 

0.030 

0.070 


fable 3.2 Beattachment 

Lengths 

Rch 

0 

»®9 

Xj^/H 

1.5 X 10“ 

0 

140 

6.5 

1.5 X 10“ 

.57 

150 

4.3 

2.2 X 10“ 

0 

180 

7.0 

2.2 X 10“ 

.57 

190 

4.5 

3.7 X 10“ 

0 

280 

6.8 

3.7 X 10“ 

.57 

25C 

5.3 
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Table 3.3 Eddy Fox«ati<m Fositioa 




(LDV) 

(Scb.) 

1.5 X 10“ 

0.0 

0.5 ♦ .3 


-.5 X 10“ 

0.57 

1.0 ♦ .3 

0.89 ♦ .14 

2.2 X 10“ 

0.0 

0.5 ♦_ .3 


2.2 X 10“ 

0.57 

l.O ♦ .3 

0.67 ♦ .08 
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Figure 3.9 Hot wire spectra of the u-velocity at the step 
without packing as used by Ganji and Sawyer 
(1980) (x/H = 0, y/H = 0.5, z/H = 0, <f = 0). 
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Figure 3.10 Hot wire spectra of the u-velocity at the 
step, with stainless steel packing and 
screens (x/H«0, y/H«0.5, z/H«0, 
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Figure 3.15 Streamwise velocity flowfield -- high Reynolds number 
(R6 u* 3.7 X lO", Uo»22.2 m/s, flame boundary). 















igure 3.18 Vorticity thickness of the non-reacttng 
shear layers. 
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Figure 3.21a Upper boundary of shear layer for 
Re^ • 3.7 X 10'. 



X/H 


Figure 3.21b Centerline position of shear layer for 
Ren*3.7 X 10'. 
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Figure 3.23 Streamwise turbulence Intensity profiles 
(Reu* 1.5 X 10\ Un= 9.12 m/s). 






119 



Figure 3.24 Streamwise turbulence intensity profiles -- medium Reyno 
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Figure 3.2G Streamwise turbuletice intensity profiles 
(Rei,= 3.7 X 10\ Uq = 22.2 m/s). 
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Figure 3.27 Transverse turbulence intensity profiles (Ren *2.2 x lO", Uo*13.3 m/s). 
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Figure 3.28 Streei»f1s€ velocity PDFs in the isotberwal layer at x/H»1.0 
(Reu» 1.5 X 10\ Uq*9.12 •/$,♦ *0). 
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Figure 3.30 Skewness factors for the low Reynolds number flow 
(Re^ * 1 .5 X 10^, Uo“9.12 m/s, — - flame boundary). 
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(Ren * 1 -5 X 10“, Uo»9.12 m/s, — - flame boundary). 
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Figure 3.36 ikewness factors for the reacting flow at higher Reynolds numbers 
(<>*0.57, — flame boundary). 
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Figure 3.37 Flatness factors fo- the reacting flow at higher Reynolds numbers 
(♦ = 0.57, — flame boundary). 
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Figure 3.38 Streamwlse velocity PDFs in the reacting layer at x/H»0.5 
(ReH»2.2xlO\ Uo» 13.3 m/s, <|»«0.57). 

















Sequence of schl ieren photographs showing 

coalescence (Ren -2.2 * 10“ . i - 0.57, 

0.3 ms apart). 




x/H 0 1 2 0 12 


figure 3.42 Sequence of schlieren photographs at the high 

entrance Reynolds rnwiber (Rch =3.7 x 10" , $ = 0.5?, 
Tgxp = 58ps, 0.6 ms apart). 
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Figure 3.46 Frequency spectrum of the v-velocity In the 
non-reacting shear layer. 
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Figure 3.47 Passing frequency of the large scale structures in the reacting shear layer. 
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gure 3.48 Passing frequency of the large scale structures in the reacting flow beyon 
reattacbment (x/H=7, (j» = 0.57). 
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Figure 3,49 


Sequence of lDV spectra takeM In 
the reacting flow (ReH* 2.2 x O'*, 
4 ■ 0,57, V - velocity componer 
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Figure 3.50 Strouhal number based on step height determined 
from schlleren and LDV measurei..ents . 
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Figure 3.52 Inverse Strouhal number based on the step 
height for the Isothermal layer. 
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Figure 3.53 Strouhal number baued on the step height 
for the Isothermal data In mixing layer 
coordinates. 
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Figure 3.54 Strouhal number based on step height for the 
reacting data in mixing layer coordinates. 
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Chapter 4 

8DIM4R AID OOlCLDSiaiS 

The priaary eia of this work is to Msess the effect of 
coabuetioQ on the rearward- feeing atep flow. The flow field ia 

docueented by laaer Doppler velociiaetry aeasureaeDts of the turbulent 
velocity and high speed achlieren visual izat ion. The results and 
conclusions of the investigation with regard to the different aspects 
of the flov field are given in the following sections. 

Initial Conditions 

The state of '..he initial boundary layer at separation is 
detensined. The boundary layer is partially tripped having high 
turbulence levels «ith a laminar profile. This partial tripping has 
t%K> effects. The linear growth rate reg >n of the shear layer is 
increased and the virtual origin appears dotmstream of the step. 

Mean Velocity Flow Field 

The mean velocity flow field is measured for both reacting 
non-reacting flow. The primary effects of combustion are: 

1. The virtual origin of the reacting shear layer defined by the 
mean velocity profiles appears farther downstream of the step but 
the growth rate of the shear layer is unchanged. 

2. In the reacting shear layer, the flame as seen from schlieren 
observations propagates faster into the incoming inviscid 
reactants than the shear layer defined by the swan velocity. 

3. The reacting shear layer defined by the mean velocity is shifted 
toward the recirculation zone. Heat expansion from the reacting 
eddies propagating into the premixed reactants increases the 
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velocity «t the top of the leyer end lowera the upper boundary of 
the cheer leyer defined by the aean velocity. Heat expanaion 
increaaea the velocity difference in the lower edge of the ahear 
layer and caucea it to grow into the recirculation none aore 
rapidly. 

4. The reattachaent length of the reacting ahear layer ia ahortened 
by 30Z due to the ahift of the ahear layer toward the 
recirculation xone. The tranaverae velocity (v) alao ahowa thia 
downward ahift. 

Turbulence Oevelopawnt 

The turbulence intenaity and the higher order aoaenta of the 

streaawiae turbulence are aeaaured throughout the flow field with the 

following reaults: 

1. Combustion causes the peak positions of the turbulence intensity 
profiles to shift toward the recirculation sone confirming the 
mean velocity results. 

2. The turbulence levels decrease rapidly 1-2 step heights before 
reattachment for both reacting and non-reacting cases due to 
curvature effects. 

3. In the reacting layer, the maximum turbulence levels are 30Z 
higher and have a double-peak in their variation with x. 

4. Koi the shear layer defined by mean velocity, the non-Gaussian 
behavior is similar for both the non-reacting and reacting 
conditions. 

5. Under combustion, there is non-Gaussian behavior above the shear 
layer defined by the mean velocity. The velocity PDF is 
positively skewed from the impelling effect of large scale 
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•tructurea «cp«nding batmen the top of the eheer layer and the 
upper flaae boundary. 

6. The atreamwiee uelocity PDFa are double>peaked in the reacting 
ahear layer near the atep. The PDFa are eingle-peaked in all 
other caaea. Intenaittent foraation of the large acale 
atructurea near the atep in the reacting layer could have cauaed 
the double-peak. 

Large Scale Structure Developaent 

The large scale atructurea are observed with high speed schlieren 

photography and frequency analysis of the LDV signal. The following 

observations are aade: 

1. As reported by Ganji and Sawyer (1980), the large scale 
structures are found to control the overall structure of the 
reacting shear layer. Reaction occurs within the eddies as they 
entrain premised reactants and hot products. The growth of the 
large scale structures is tied to the propagation of the flame. 

2. The formation of the eddies is delayed in the reacting shear 
layer as the local Reynolds number is reduced by combustion. 

?. As ceen from the schlieren visualization, the propagation of the 
flawc into the incoming reactants extends above the shear layer 
defined by the mean velocity profiles. The upper edge of the 
reacting large scale structures defines the flame boundary. 

4. The most striking effect of combustion is the reduction of eddy 
coalescence in the reacting shear layer. The pairing process is 
nearly eliminated as a growth mechanism in the reacting layer. 
The layer grows primarily by entrainment of fluid and heat 
expansion. 
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